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ABSTRACT 

We explore extensions to the standard ACDM cosmological model using new measurements of the 

cosmic microwave background (CMB) from the South Pole Telescope (SPT). Adding SPT measure- 
ments to WMAP7 significantly improves constraints on possible extensions to the ACDM model; the 
addition of low-redshift measurements of the Hubble constant (Hq) and baryon acoustic oscillation 
feature (BAO) leads to further improvements. Before combining these datasets, we check for con- 
sistency in the ACDM model between measurements of the CMB (SPT+ IF^AMP7), Hq and BAO, 
and find evidence for some tension between the datasets. Within the CMB data alone, we find only 
weak support for physics beyond the ACDM model due to a slight trend of decreasing power at 
smaller angular scales, relative to the prediction of the ACDM model. This trend could be due to a 
logarithmic scale dependence of the power-law index of the primordial power spectrum, dns/dlnk. 
Alternatively, the trend coiild arise either from adjustments at small or large scales. The power at 
small scales is sensitive to the damping scale which is influenced by both the helium abundance, Yp 
and the effective number of neutrino species, A^eff- The power at large scales is affected by the inte- 
grated Sachs- Wolfe effect which is sensitive to the sum of neutrino masses, '^m.,^. These extensions 
have similar observational consequences and are partially degenerate when considered simultaneously. 
These degeneracies can weaken or enhance the apparent deviation of any single extension from the 
ACDM model. Of the 6 one-parameter model extensions considered, we find the CMB data to have 
the largest statistical preference for running with —0.046 < dng/dlnk < —0.003 at 95% confidence. 
This preference for dns/dlnk < strengthens to 2.7a for the combination of CMB+BAO+ifo- Run- 
ning of this magnitude is difficult to explain in the context of single-field, slow-roll inflation models. 
When varying the effective number of massless neutrino species, wc find N^g = 3.62 ± 0.48 for the 
CMB data. Adding Hq and BAO measurements tightens the constraint to A'off = 3.71 ± 0.35, 1.9 cr 
above the expected value for three neutrino species. Larger values of Ncs relieve the mild tension 
between the CMB, Hq, and BAO datasets in ACDM. Instead fixing A^eff to the standard value and 
allowing the sum of neutrino masses to be free, we find the combined data from the CMB, BAO, Hq, 
and Sunyaev-Zel'dovich selected galaxy cluster abundances prefer nonzero neutrino masses at 3.0 cr. 
The median value is = (0.32 ± 0.11) eV, a factor of six above the lower mass bound set by 

neutrino oscillation observations. All datasets except Hq show some preference for massive neutrinos; 
data combinations including Hq favor nonzero masses only if BAO data are also included. We also 
constrain the two-parameter extensions N^s + rriv and A^eff -l- to explore constraints on sterile 
neutrinos and big bang nucleosynthesis respectively. 

Subject headings: cosmic microwave background, cosmological parameters, early Universe 
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model conditioned on CMB measurements at I < 1000 



1. INTRODUCTION 

Measurements of the damping tail of the cosmic mi- 
crowave background (CMB) power spectrum provide new 
insights into the spectrum of the primordial density fluc- 
tuations, the contents of the pre-recombination plasma, 
and, through the effect of gravitational lensing, the low- 
redshift properties of the universe. In this paper, we 
study the constraints that can be placed on extensions to 
the standard ACDM cosmological model by recent South 
Pole Telescope (SPT) CMB power sp ectrum measure- 
ments presented in a companion paper ( Story et al.|2012" 



hereafter S12) and complementary lower-redshitt probes 
For each extension, we discuss both the cosmological con- 
straints and the physical origin of these constraints. 

As shown by S12, the SPT bandpowers lead to con- 
straints on the six standard ACDM parameters that are 
consistent with and comparable to thos e from WMAP7 



(Larson et al. 2011 
is especially remar 



Komatsu et al. 2011 ). The agreement 
table given that the smaller angular 



scales probed by the SPT data are sensitive to new phys- 
ical effects that are unimportant on the larger angular 
scales probed by WMAP. The combination of SPT and 
WMAP7 data can be used to exploit these additional ef- 
fects in order to improve constraints on the ACDM model 
and possible extensions. 

Measurements of the damping tail of the CMB power 
spectrum are sensitive to the low-redshift universe 
through the gravitational lensing of the CMB as well 
as the angular-diameter distance to the last scattering 
surface. The lensing sensitivity introduces a second 
constraint on the late-tim e expansion rate that breaks 



the geometr ic degeneracy ( Bond et al. [I997 Zaldarriaga 
et al. 19971 that exists in large-scale CMB data. Mea- 



surements of CMB lensing, like that recently reported by 
S12, can be used to constrain the mean curvature of the 
observable universe. 

More importantly, the SPT data expand the angu- 
lar range over which the CMB temperature power spec- 
trum is probed, leading to improved measurem ents o f the 
damping scale due to photon diffusion (Silk 1968) and 
the positions of acoustic peaks. Past CMB experiments 
have found evidence for a trend of decreasing power at 
high multipoles, relative to the predictions of the ACDM 
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( Harnann et al.|2010b Dunkley et al.|20TT Keisler et al 
201 1| ), a trend that persists with the inclusion ot the tuU- 



survey SPT bandpowers. This trend may be thought 
of as a scale-dependent tilt in the CMB power spec- 
trum, and we explore five one-parameter extensions to 
ACDM that effectively allow for such a scale-dependent 
tilt. These extensions are: allowing massive neutrinos 
(^TOi,), introducing "running" of the spectral index of 
the primordial power spectrum (drtj/dln fc), varying the 
effective number of neutrino species (Ncg), allowing the 
helium abundance (Yp) to depart from the predictions of 
big bang nucleosynthesis (BBN), and allowing non-zero 
early dark energy (ile)- A unifying theme in this paper 
will be a preference for model extensions that can accom- 
modate a scale-dependent tilt that becomes increasingly 
red at higher multipoles. 

The CMB damping tail measurements do not add sen- 
sitivity to the dark energy equation of state w, since dark 
energy becomes important at late times and this param- 
eter has little effect on the lensing amplitude. Therefore, 
we do not present constraints on ACDM-l-u; in this pa- 
per, although we sometimes include w when exploring 
the effect of parameter degeneracies on the constraints 
for other parameters. 

We explore two 2-parameter extensions to the ACDM 
model which are physically and theoretically motivated: 
ACDM-|-7Veff+>"p and ACDM+iVeff +^ m^. The first 
case is an interesting consistency test of BBN. Both Yp 
and iVofi' are primarily constrained by measurements of 
the damping scale, however both can be constrained si- 
multaneously since varying N^s has the secondary effect 
of inducing a phase shift in the acoustic oscillations. The 
second extension, ACDM+Ncs+J^^i^i is a test of the 
possible existence of sterile neutrinos, which have been 
postulated with masses in the eV range in orde r to ex- 
plain some neutrino oscillation anomalies 



Aguilar 

et al.||200l| [Aguilar-Arevalo et aL]|2010[ [Mention eFaT 

miT. 

hi 



et al. (2012|) BAO measurement at z = 0.57 and the 



While all of these extensions can be constrained with 
CMB data alone, low-redshift measurements add con- 
straining power and provide consistency checks to the 
CMB data. We incorporate measurements of the Hubble 
constant (Hq), the baryon acoustic oscillation (BAO) fea- 
ture of the matter power spectrum, luminosity distances 
from supernovae (SNe), the matter power spectrum in- 
ferred from a redshift catalog of luminous red galaxies 
(LRGs), and the abundance of Sunyaev-Zel'dovich (SZ)- 
selected clusters. We pay particular attenti on to the 
BAO and H p measurements, specifically the Anderson| 



Kiess et al. "( [2011 ) determination of Hq, because they are 
highly precise constraints that interact in an interesting 
way with those from the CMB. We quantify the level of 
consistency between the CMB, BAO, and Hq datasets in 
the context of the ACDM model, and we study which 
model extensions are preferred by the CMB-only and 
combined CMB+BAO-f^iJo datasets. 

Throughout this work, we make a point of explain- 
ing the physical mechanisms behind the constraints we 
find. That is, we address questions such as, "-ffow; are 
the neutrino masses constrained?" and " What features 
in the data explain the preference for nonzero masses?" 
We also clarify what other extensions might explain the 
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same features in the data. 

The structure of the paper is as follows. We briefly de- 
scribe the analysis methodology and data from the SPT 
and other experiments used in this analysis in § [2] The 
physical origins of the CMB constraints on the standard 
ACDM cosmology are presented in § [3j We discuss the 
consistency between datasets, and the evidence for exten- 
sions to the ACDM model in § [4] We study the physical 
origins of the constraint on non-zero curvature in § [51 
We consider massive neutrinos in § [6] and running of the 
scalar index in § [7j Other single-parameter extensions 
- iVoff , Yp , and ~ are discussed in § [8j We consider 
2 -pa rameter extensions in § [9j Finally, we conclude in 



2. PRELIMINARIES 

Here we introduce the datasets that will be used in 
the cosmological parameter fitting. We also describe the 
Markov Chain Monte Carlo (MCMC) methods that will 
be used to fit cosmological parameters. 

2.1. CMB Power Spectrum Measurements 

In this work, we will use CMB bandpowers from 
WMAP and SPT, and will refer to the combined dataset 
as the "CMB" dataset. The WMAP data we use are 
the measurements of the CMB power spe ctrum at 2 < 
I < 1 200 from the WMAP 7-year release (Larson et al. 
20111.^® We will focus on the new SPT power spec- 



trum results presented by S12, hereafter referred to as 
the SPT dataset. The SPT bandpowers represent a two- 
fold improvement o ver the best previous measurement 
(Keisler et al. 2011 Kll) of the CMB power spectrum 
at / > 800. We briefly summarize the analysis of S12 
and review the modeling of the contributions from ex- 
tragalactic foregrounds. 

The full SPT-SZ survey was completed in November 
2011 and covers ^2500 deg^ in three frequency bands at 
95, 150, and 220 GHz. The SPT bandpowers used in this 
work are based on the 150 GHz data from the complete 
survey. As in earlier SPT power spectrum analyses, the 
bandpowers are calculated using a pseudo-C; method. 
The bandpowers cover the range 650 < I < 3000 and are 
sample-variance limited at I < 2900. The covariance ma- 
trix contains terms for the sample and noise variances as 
well as beam and calibration uncertainties. Beams are es- 
timated through a combination of observations of planets 
and bright synchrotron sources in the survey fields. The 
bandpowers are calibrated with an accuracy of 2.6% in 
power by comparing the SPT and WMAP 7-year power 
over the multipole range I £ [650, 1000]. 

We model contributions from extragalactic foregrounds 
to the SPT bandpowers exactly as in S12. Using the 
notation Di = l{l + l)C;/(27r) these are given by 
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7-)SZ rSZ 
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(1) 

where D^q^q, D^^q, and i'fooo are the amplitudes at I — 
3000 for the Poisson point source power, clustered point 

^® A few days prior to su bmi ssion, the WMAP collaboration 
released the 9-year data I jHinshaw et al. 2012 Bennett et al. 2012fc. 
Due to time constraints a!nd the consistency of the y-year data with 
the 7-year data, we chose not to incorporate the new data into this 
work. 



source power, and total SZ power spectrum respectively, 
and /f ^ gives t he ^-dependence of t he SZ power spectrum 
as modeled in Shaw et al. (2010|). We apply Gaussian 
priors on these amplitudes of D^qqq — 19.3 ± 3.5 nK^ 



j :)CL„ = 5.0±2.5A <K^ and ff^go = 5.5 ± 3.0 /xK^ (|Shi- 



rokoff et al. 



^3000 , , 

2011[ ). We include the SZ term inTlie 

WMAPl likelihood calculation, but we do not include 
the two point source terms since the WMAPl likelihood 
code already marginalizes over appropriate point source 
terms. We have tested that all parameter constraints are 
insensitive to loosening the assumed extragalactic fore- 
ground priors. 

As was argued by S12, we expect galactic foregrounds 
to have negligible impact on the SPT bandpowers and do 
not include a contribution from them in the band power 
modeling. 

2.2. External Datasets 

Low-redshift observations can substantially inform the 
constraints on some cosmological parameters, particu- 
larly those parameters with impact at late times such as 
the mean curvature of the observable universe. We in- 
clude these low-redshift datasets in many of the cosmo- 
logical constraints presented in this work. These datasets 
are described below. 

2.2.1. Baryon Acoustic Oscillations 

BAO measurements use the spatial correlation of 
galaxies to constrain the comoving size of the sound 
horizon at the epoch when baryons and photons de- 
couple, rs(zdrag)- This redshift, Zd^ag: is defined by 
r&(77(zdrag)) = 1, where 



3p6 



(2) 



Here ry is conformal time with 770 the conformal time to- 
day, a is the scale factor, pf, and are the baryon and 
photon densities, is the number density of free elec- 
trons, and ar is the Thomson cross section. The BAO 
feature in the galaxy correlation function can be observed 
as a function of both angular and redshift separations. 
Thus galaxy surveys can, in principle, constrain both 
DA{z)/rs{zdYe,^) and (z)rs(zdrag)- However, the cur- 
rent generation of BAO observations typically reports 
constraints on a hybrid quantity, T's(zdrag)/^y(2cff) 
(henceforth rg/Dy), where Dy{z) — D\{z)cz / H {z) ^ c 
is the speed of light and ZeS is the effective redshift for 
the population of galaxies in the survey. 

We use measurements of rs/Dy from the followi ng ex- 
perim ents: SDSS at Zcff = 0.35 ([Pa dmanabhan et al. 
20T2I , BOSS at Zeff = 0.57 ( [Anderso n ct al. 2012)7and 
WTggleZ at z^g = 0.44,0.6, and 0.73 {Blake et ai.,|20TT| ). 
All of these experiments report rg/Dy (or its mverse) us- 
ing th e fittin g formula for Zdrag presented by Eisenstein 
fc Hu[ (1998). We numerically calculate ^^(zdrag j with 
CAMBl [Lewis et al.l2000[ ). To account for the difference 
between Iculated by the fitting formula with 



•Zdrag as calculated by CAMB, Hamann et al. (2010a I 



find it necessary to r escale rg by a factor 154.66/150.82 
Mehta et al.| (20121 show this constant rescaling ade 



quately accounts tor the differences in Zdrag over the 
range of baryon density and matter density allowed by 
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WMAP7. We therefore apply this rescaUng to the re- 
ported Ts/Dv data for all parameter fitting. 

2.2.2. The Hubble Constant 

In contrast to the BAO approach to measuring cos- 
mological distances, which starts from a known distance 
scale at z ^ 1100, measurements of the Hubble con- 
stant are based on the classic distance ladder which is 



built outwards from the solar system. Riess et al. (2011 1 
use 253 Type la SNe to reach out into the Hubble flow 
and calibrate their luminosity distances with over 600 
Cepheid variable stars in the host galaxies of eight nearby 
Type la SNe. The Cepheids are calibrated in turn with 
three different methods, including 13 with trigonometric 
parallax. Accounting f or both statistical and systematic 
sources of uncertainty, Riess et al. (20111 conclude that 



Hp = 73.8 ± 2.4 km/s/M pc 



Freedman et al. (2012) recently released a new analysis 
01 the HST Kc5rProject, with a Spitzer-ha.sed calibration 
of the Cepheid distance scale, in which they find Hq — 
74.3±2.1 km/s/Mpc. As we had already completed many 
of the Markov chains used in this work, we chose not to 
rerun with this slightly improved constraint. We note 
that the tension in the ACDM model between BAO and 
Hq discussed later would mildly increase with this new 
data since the measured Hq value increases by 0.2 a, and 
the uncertainty decreases by ^ 10%. 

2.2.3. Luminosity Distances from Supernovae 

Type la SNe, employed as standardizable candles, re- 
main the most precise means to extend the distance lad- 
der to redshifts that are sensitive to cosmic acceleration. 
We us e a sample of 427 SNela assembled by |Conley et al.| 
(2011) to constrai n the luminosit y dista nced We follow 
the treatment in Sullivan et al. (2010) to handle the 
dependence of supernovae absolute magnitude on host 
galaxy stella r mass. We follow the approach of |Conley| 
et al. (2011) to marginalize over systematic uncertain- 
ties related to the uncertain intrinsic properties of su- 
pernovae. 

2.2.4. Matter Power Spectrum from SDSS DR7 LRG 



Reid et al. (2010) use a catalog of luminous red galax- 
ies trom the seventh data release of the Sloan Digital Sky 
Survey (SDSS DR7) to reconstruct the dark matter halo 
density field at z < 0.4 over nearly 8,000 deg^. They use 
the reconstructed density field to determine the matter 
power spectrum over the range 0.02 < k < 0.2 /iMpc"^. 
We use their published likelihood code to incorporate 
this information in cosmological constraints. The in- 
formation is primarily from the shape of the spectrum 
since the unknown halo masses, and therefore unknown 
biases, prevent the amplitude of the measured matter 
power spectrum from being used. 

2.2.5. SPT SZ-selected Clusters 



Reichardt et al. (2012b R12) present a catalog of 224 
galaxy cluster candidates selected by their SZ signifi- 
cance from the first 720 deg^ of sky surveyed by the 
SPT. As in R12, we use only the 100 cluster candidates 
at z > 0.3 with a significance greater than five. The 
optical properties of these systems are detailed by |Song 



et al. (2012). X-ray Yx measurements of 14 of the clus- 
ters are also used in the cosmological constraints; these 
measu reme nts have been described by [Andersson et al. 
(|2011| andlBenson et al.l (|2011|). 



2.3. Parameter Estimation Methodology 

Posterior probability distributions for the parame- 
ters are calculated using Markov Chain Mo nte Carlo 
(MCMC) methods fChr istensen et al.||200T|) with the 
publi cly available CosmoMC^' package (Lewis & Bridle 
2002|). We have mod ified the code to include the fore- 
2.11 and the extensions to ACDM 



ground terms from 
detailed in later sections; these modules and instructions 
for compiling them are available at the SPT website.^* 
We use the high_accuracy_default option in CosmoMC, 
and by default call the January 2012 version of CAMB 
(Lewis et al. 2000) to calculate the CMB and matter 
power spectral Unless otherwise noted (for instance in 



6.2 ) , we will report parameter constraints based on the 
median value of the posterior distribution with 1 a error 
bars defined based on the interval that includes the cen- 
tral 15.85% to 84.15% of the MCMC samples. For cases 
when we compare a parameter to some fiducial value 
(e.g. an extension parameter to its ACDM-consistent 
value), we report-unless otherwise noted-a significance 
in a based on the Gaussian z-score of the probability 
above (or below) the fiducial value. 

We will present parameter constraints on a spatially 
flat, gravitationally lensed ACDM model along with se- 
lected one- and two-parameter extensions to this model. 
The properties of the standard ACDM model are spec- 
ified by six parameters: the optical depth to reioniza- 
tion (r), the baryon density (fibft,^), the cold dark mat- 
ter density (fic/i^), the angular scale of the sound hori- 
zon at last scattering (Og), the scalar spectral index (n^), 
and the amplitude of primordial scalar perturbations at 
wavenumber k — 0.05 Mpc~^ (^/j)- Throughout this 
paper, we will use the notation ujx = fl^h^. Since 
Qx = Px/Pc where pc = 3Hq / (8ttG) is the critical den- 
sity, ujx is independent of the Hubble constant. We define 
the total matter density as the sum of baryons and cold 
dark matter, so that ujm = ijJb + ^c- AH densities are 
evaluated at the present epoch. 

To speed up parameter estimation, we make use of 
the Boltzma nn code emulator PICO (Fendt & Wan- 
delt||2007b|a[ ). PICO works by empirically modeling the 
CMB and rnatter power spectra as sums of polynomi- 
als, with the (cosmological parameter-dependent) poly- 
nomial coefficients determined by fitting them to output 
from CAMB. We have trained PICO for a 10-parameter 
model that includes the six ACDM parameters as well 
as Yp, Ncs, ^m^, and dus/dlnk. We use PICO in- 
stead of CAMB when working with any subset of this 
model space. Evaluating the power spectrum with PICO 
is several orders of magnitude faster than doing so with 
CAMB, and is even more accurate due to its training 
on CAMB runs with higher numerical accuracy than the 
high_accuracy_default setting. The PICO code and re- 
sults of this training are available at the PICO website. 

http://cosmologist.info/cosmomc - January 2012 version 
littp:/ /pole. uchicago.edu/public/data/storyl2/ 
https : //sites . google . com/a/ucdavis . edu/pico] 
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Fig. 1. — This figure illustrates the agreement between parameter 
values derived from the SPT and WMAP7 bandpowers. For each 
ACDM parameter, we plot the fractional change in the parameter, 
defined by (asPT/c«WAfAP7 — l)i and the error bar for this quan- 
tity. Following S12, the amplitude of scalar perturbations, A|j, 
is evaluated at fc = 0.05 Mpc~^. The optical depth to reioniza- 
tion is not plotted, as a WMAP7-ha.sed prior is used in the fits to 
the SPT bandpowers. The most significant shift is to the sound 
horizon scale which increases by 1.5 cr. 



3. ACDM RESULTS 

In this section, we consider the SPT-only and 
WMAP7-on\y constraints on the ACDM model param- 
eters. As shown by S12 and Figure [T] the parameters 
estimated from the two datasets are consistent. Here we 
focus on the physical origins of the baryon and matter 
density constraints derived from different angular scales 
probed by WMAP7 and SPT. For each case, we begin by 
briefly re viewing the origins of the WMAP7 constraints 
(see, e.g., Hu fc Dodelson|2002 ) and then discuss the ori- 
gins of the yP i' constraiiits. As we will see in later sec- 
tions, understanding these mechanisms provides a foun- 
dation from which to explore other datasets and model 
extensions. 

The WMAP7 constraint on the matter density uim de- 
rives from differences in the acoustic oscillations during 
the radiation-dominated and matter-dominated epochs. 
Modes that "enter the horizon" "^^ during radiation domi- 
nation receive a boost to their amplitude from the decay 
of the gravitational potential that occurs during their 
first compression into potential wells. In an expanding 
universe, the plasma would need to collapse at the free- 
fall (pressureless) rate to preserve the potential. The 
pressure support of the plasma slows the collapse, leading 
to a decaying gravitational potential. The decay boosts 
the oscillation amplitude since modes gain more energy 
during the compression into the initially larger poten- 
tial well than they later lose during the expansion out 
of the much reduced potential well. Modes entering dur- 
ing matter domination do not see such a boost, since 



the potential is predominantly sourced by pressureless 
dark matter. The difference in the oscillation ampli- 
tudes allows us to determine the angular scale separating 
modes that entered the horizon during matter domina- 
tion from those that entered during radiation domina- 
tion, ^oqj which can be related to the redshift of matter- 
radiation equality, Zgq. Because the radiation density, 
Wj-ad, is completely specified in the ACDM model by the 
well-known temperature of the microwave background 
today and assumptions about neutrino production, this 
measurement of z^q corresponds to a determination of 
uj„i through the relation 1 -I- z^q = Wm/wrad- 

The angular scales probed by the SPT bandpow- 
ers correspond to modes that entered the horizon in 
the radiation-dominated era, since the transition occurs 
around Icq — ir/dcq — 434. Thus the SPT bandpowers 
are largely insensitive to ajj^ through the mechanism de- 
scribed above. To confirm this, we have tested removing 
the largest scales (third peak, I < 1000) from the SPT 
bandpowers and find only a minimal (5%) degradation 
of the SPT-only ojm constraint. The SPT data are in- 
stead sensitive to the amplitude of gravitational lensing 
which in turn depends on aj„j. Gravitational lensing of 
the CMB smooths out the peak-trough structure of the 
CMB power spectrum, allowing lensing to be detected in 
the SPT power spectrum. The magnitude of the lensing 
potential power spectrum increases with uJm- To sup- 
port the claim that the SPT bandpowers determine Wm 
through its effect on the amplitude of the lensing power 
spectrum, we have tested removing the lensing amplitude 
information. We accomplish this by rescaling the lensing 
power spectrum, Cf"^, by a free parameter, A^, defined 
by: 



(3) 



The horizon expands with time; a Fourier mode enters the 
horizon when its wavelength equals the horizon size. 



Marginalizing over Aj^ doubles the uncertainty on w^, 
degrading the constraint from LUm = 0.1286 ± 0.0071 to 
UJm = 0.129 ± 0.013. When SPT and WMAP7 data are 
combined, uJm is primarily constrained through its im- 
pact on ^cq- 

The WMAP7 constraint on Wf, derives from the relative 
heights of the even and odd acoustic peaks. The baryon 
density affects the plasma's pressure, which impacts the 
acoustic oscillations. Changing the pressure has opposite 
effects on compression into gravitational potential wells 
(odd peaks) versus compression onto gravitational poten- 
tial hills (even peaks). The net result is that increasing 
uJh raises the odd peaks relative to the even peaks. 

The SPT constraint on baryon density is primarily de- 
rived from the same even/odd peak modulation. Due to 
the decay of the gravitational potentials, the even/odd 
modulation of peak heights is suppressed on smaller 
scales. The dominant constraining power in the SPT 
data is from the height of the third acoustic peak. In- 
creasing Uf, at fixed Zoq and 0s increases the height of the 
third acoustic peak relative to higher peaks. If we re- 
move the SPT bandpowers over the third acoustic peak 
(dropping all bandpowers at I < 1000), the SPT-only 
constraint on cjf, degrades by nearly a factor of 2 from 
ojb = 0.0230 ± 0.0011 to ojb = 0.0236 ± 0.0020. 

The damping scale provides a second (weaker) lever 
arm to constrain the baryon density from the SPT 
data. Photons are imperfectly coupled to baryons in 
the primordial plasma, and the resulting photon diffu- 
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sion damps the acoustic oscillations. This photon dif- 
fusion is described by the diffusion length r^, the root 
mean squared comoving distance a photon has traveled 
over the time in which the scale factor grows from to 
a* (the scale factor at decoupling). The diffusion length 
is given approximately by 



da 



(4) 



Here is the number density of free electrons, ax is the 
Thompson cross-section, and R = 3pb/(4p^). The factor 
in square brackets is due to the directional and polariza- 



tion dependence of Thomp son scattering (Kaiser 1983 
Zaldarriaga fc Harari||1995[ ). The diffusion length is af- 
fected by uJb through the number density of free electrons, 
nf,(z) cx Xe{z)u}h (1 — Yp), where Xe{z) is the fraction of 
Hydrogen atoms that are ionized and Yp is the primordial 
fraction of baryonic mass in helium. 

We quantify the importance of the SPT damping scale 
measurement in constraining the baryon density by in- 
troducing Yp as a free parameter. At fixed Wf,, vary- 
ing Yp only impacts the free electron density and hence 
the diffusion scale. When we allow Yp to vary freely, 
the uncertainty on Wfc slightly increases from 0.0011 to 
0.0013. Thus, the primary sensitivity of the SPT data to 
the baryon density is from the measurement of the third 
acoustic peak rather than damping scale. 

The SPT constraints on the other ACDM model pa- 
rameters are based on similar physical mechanisms to 
the WMAP7 constraints. The angular size of the sound 
horizon, 6^ is constrained via the peak locations, while 
the primordial parameters Ug and A|j^ are constrained 
through their impact on the amplitude and tilt of the 
CMB power spectrum. As mentioned in S12, the SPT 
bandpowers do not constrain r; we instead apply a prior 
based on the WMAP7 polarization measurements. 

We note that there is a trend in the CMB data for de- 
creasing power at higher multipoles, relative to fits to the 
ACDM model using data from I < 1000. In the ACDM 
model space, this deficit manifests as the downward shift 
in Us between WMAP7 and SPT seen in Figure [ij As 
we shall see in later sections, model extensions that can 
accommodate this trend are somewhat preferred by the 
CMB data. 

4. CONSISTENCY OF CMB, BAO, AND Hq 
MEASUREMENTS AND EVIDENCE FOR EXTENSIONS 
TO THE STANDARD COSMOLOGICAL MODEL 

Combining the BAO and Hq measurements with the 
CMB will have important implications for the model ex- 
tensions we consider. Here we investigate the consistency 
of these datasets with one another. Tension between 
two sets of data in the context of a ACDM model can 
be caused by (1) statistical fluctuations in the measured 
values, (2) an inaccurate or incomplete accounting of the 
uncertainties (or biases) in one or more of the datasets, 
or (3) the inadequacy of the cosmological model. We dis- 
cuss the consistency between datasets in § |4.1| and discuss 
any evidence for extensions to the ACDM model in § |4.2[ 

4.1. Consistency between datasets 

In the context of the ACDM model, we quantify the 
consistency between two datasets, x and y, with the 



TABLE 1 

Relative Consistency Between Datasets in the 
ACDM Model 



Dataset Comparison 


Add. 


Ax' 


corrected 




dof 




PTE 


(CMB) vs i+Ho) 


1 


0.23 




(CMB) vs (+BAOsDSs) 


1 


1.8 




(CMB) vs (■+BAOwiggioz) 


3 


1.7 




(CMB) vs (+BAOboss) 


1 


5.1 


0.090 


(CMB-t-i^o) vs (-I-BAO) 


3 


7.4 




(CMB-I-BAO) vs i+Ho) 


1 


3.5 


0.12 



Top 4 rows: the relative consistency between the 
CMB and the individual low-redshift measurements. 
For each model, we calculate the Ax' between the 
x' ■ for the CMB dataset and the x' ■ for the 
CMB-|-(each low-redshift experiment). The least con- 
sistent, BAObosSi is shown in row four, along with 
the probability to get a larger Ax' given four dataset 
combinations, specifically: (corrected PTE) = 1 — 
(1 - PTE)'^. 

Bottom 2 rows: the Ax' between the CMB-t-i^o 
and CMB-|-BAO+-ffo datasets combinations, and the 
CMB-I-BAO and CMB-|-BAO-|-//o combinations. For 
the least consistent comparison, (CMB-fBAO) vs 
{+Ho), we show the probability to get a larger Ax' 
given two chances. 



use of a symmetric Ax^ statistic. We calculate Ax^ = 
xl+y — xl — xl where xl+y is the minimum value of the 
X^ obtained for the joint (x + y) dataset in the model M 
with a similar definition for x^ s-nd x^. We use a mod- 



ified simulated annealing minimizer (Goffe et al. 19941 
to find the minimum x^ point. The Ax'^ is proportional 
to a likelihood ratio statistic, and is also connected to a 
Bayesian evidence ratio in t he case of normal prob ability 
distributions, as shown by Marshall et al. ( |2006 ). The 



expected distribution of Ax^ is a x'^ 



distribution with 



degrees of freedom (dof) N^^2 = N^+y — — Ny. Here 
Nx is the dof for dataset x. We then compare the ob- 
served value of Ax^ with its expected distribution and 
find the probability to exceed (PTE) this observed value. 

For example, to compare the consistency of the 
CMB dataset with the Hq "dataset," we calculate 

Xmin,[CMB-|-Ho] ~ ^min,[CMB] ~ Xmin,[Ho] BcCaUSC 

Hq has only one datapoint while the model has six 
free parameters, we have Nng = 0, Ncmb = X, and 
NcMB+Ho = X + 1. Thus Ax^ in this comparison should 
be drawn from a distribution with 1 dof, and we find 
a PTE of 0.63. 

As we are considering the consistency between the 
CMB and multiple datasets, we have a set of size N 
of these PTEs, where N corresponds to the number of 
dataset comparisons in the test. The minimum PTE 
should yield the strongest evidence of possible tension. 
We therefore compare the observed minimum PTE to 
the expected distribution of the minimum of a set of N 
uniformly distributed random variables to derive a "cor- 
rected PTE". Specifically, the corrected PTE a is de- 
rived from the uncorrected PTE ao by a = 1 — (1 — ao)''^. 
The corrected PTE for Ax^ represents a measure of 
tension between datasets under the assumption that 
the given model is correct, and the difference between 
datasets is due to random scatter in the measurements. 
A high PTE indicates the data are remarkably consistent. 
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Fig. 2. — Here we illustrate the degree of consistency between the CMB data and the two external datasets which are in the most 
tension: BOSS and Hq. The solid lines mark the 1 and 2cr contours from the CMB data in the Hq — Vs / Dy (0-57) plane, while the color 
encodes the value of Clmh^, as shown in the color scale on the right. From left to right, the CMB data used is WMAP7, SPT+t, and 
SPT+WMj4P7. The horizontal solid and dashed lines mark the central value and la region for the BOSS BAO measurement, while the 
vertical lines do the same for Hq. The joint 1 and 2 a likelihood regions for BOSS+//o measurements are denoted by the dark and light 
grey shaded contours. 



and low PTE indicates the presence of tension between 
the measurements. We note however that if the model is 
poor, the corrected PTE may not be a good representa- 
tion of tension with the model. 

In the top four rows of Table [ij we report the Ax^ be- 
tween the CMB and each individual low-redshift dataset 
in the ACDM model. Among the BAO datasets, the 
BAOwiggicZ dataset has a low Ax^ of 1.7 (3 dof), while 
the BAOboss dataset has an unusually high Ax^ of 5.1 
(1 dof), driven in part by the high precision of this data 
point. The probability of getting a Ax^ this large given 1 
dof (un-corrected PTE) is 2.3%; in the case of 4 random 
draws (corresponding to the 3 BAO plus 1 Hq datasets) 
from this Ax^ distribution the probability of at least one 
being this high (corrected PTE) is 9.0%. 

When the three BAO measurements are combined into 
a single dataset, the Ax^ contributions from BAOwiggicZ 
and BAOboss balance out to give a more reasonable to- 
tal Ax^. Thus the CMB-I-BAO versus Hq case has a 
lower PTE than the CMB+Hq versus BAO case. The 
Ax^ between CMB-hBAO and Hq is 3.5 for 1 dof, and the 
corresponding probability of getting a lower PTE given 
two random draws is 12%. 

From Table m we see some degree of tension between 
the CMB and the BOSS BAO data point, and between 
CMB -I- BAO and Hq. Both of these tensions exist (to 
a lesser degree) with the WMAP7 data alone and have 
been no ted by others (Anderson et al.|2012 Mehta et al. 
2012| . [Anderson et ai.| ^2U12"| highlighted the impor- 
tance or ujm for consistency between WMAPl and BAO 
measurements, noting that a la increase in io„i from 
its WMAPT-'mieiied value would bring them into near- 
perfect agreement. Instead, the SPT-I- WMAP7 data pre- 
fer a slight downward shift in ojm which, together with 
the decreased uncertainty, leads to the increased level of 
tension between BAO and the CMB in the ACDM model. 
At the same time, the shift in Wm improves agreement of 
the CMB data with the Hq measurement. 

The relationship between and the CMB, BAO, and 
Hq measurements can be seen in Figure [2] This figure 



shows the parameter space Hq vs. rs/Dv{z = 0.57) (the 
characteristic BAO parameter at the redshift reported 
by BOSS), which is the plane in parameter space that 
is constrained by the Hq and BAOboss datasets. When 
combined with any of the three combinations of SPT and 
WMAP7 data shown in this figure, the Hq data prefers 
lower values of oJm, while the BAO data prefers higher 

We find the apparent tension significant enough in 
some model spaces, including ACDM, to suggest caution 
in interpretation of the results. However, in no model 
spaces is the significance sufficient to rule out statistical 
fluctuations, and we have no evidence for either system- 
atic biases or underestimated uncertainties in the data. 

4.2. Evidence for Extensions 

Two natural questions are: (1) which extensions to 
the ACDM model are preferred by the data, and (2) do 
these extensions ease or increase the tension between the 
datasets? Here we answer these two questions, then give 
physical explanations for the results. 

We define a slightly different Ax^ statistic to quantify 
how strongly the data prefer a given extension (f> over 
the ACDM model: Ax^ = Xacdm " Xacdivh-0- A large 
reduction in Xmin (equivalently a large Ax^) relative to 
the baseline ACDM model indicates that the extended 
model is a much better fit to the data. For instance, 
if Ax^ = 4 (6.2) for a given extension with one (two) 
degrees of freedom, this means that the data favor this 
extension at 2cr. We consider models which are favored 
by more than 2 tr to be "preferred" by the data. However, 
it is worth noting that we do not consider all possible 
physically motivated extensions, and also that a single 
feature in the data may lead to an apparent preference for 
multiple extensions with degenerate effects. The results 
of this statistic for the extensions we consider in this work 
are presented in Table [2j 

The CMB dataset shows a > 2 ct preference for the ex- 
tensions AGBM+dns/dlnk and ACDM-I-Fp. The Ax^ 
for ACDM-|-dn,/dlnfc and ACDM-t-Fp is 4.9 and 4.4 re- 
spectively. Allowing running or varying the helium abun- 
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Fig. 3. — This figure expands Figure [2] to investigate the consistency between the CMB, BAObosSi ^iid Hq datasets in one-parameter 
extensions to ACDM. Each panel shows the Ho — rs/Dv(0.57) plane for a different extension. The colored dots are samples drawn from 
the SPT-I- WMAP7 MCMC chain, with the color coding reflecting the value of the extension parameter, as shown in the color scales on the 
right. The horizontal solid and dashed lines mark the central value and 1 a region for the BAOboss measurement, while the vertical lines 
do the same for Hq. The joint 1 and 2cr likelihood region for BAOboss+^^0 measurements is denoted by the dark and light grey shaded 
contours. The blue (red) contours overlaid show the 68% and 95% confidence regions for CMB+BAO (CMS+Hq)- Varying the effective 
number of neutrino species leads to the best agreement between the CMB, BAObosSi ^i^d ^0- 



9 



dance allows for an increasingly red, scale-dependent tilt, 
and provides a better match to the observed CMB band- 
powers than does the ACDM model. All other single- 
parameter extensions are not significantly favored by the 
CMB data. The combined CMB+BAO+i7o data show 
even stronger preference for these same model extensions, 
with Ax^ for ACBM+dn./dlnk and ACDM+lp of 7.4 
and 5.3, respectively. 

The CMB alone prefers neither two-parameter ex- 
tension at >2cr. Although the Ax^ = 5.5 for the 
ACDM+ Nctf+Yp model is higher than that of any single 
parameter extension given the CMB data, this is for two, 
instead of one, additional dof. The ACDM-t-A^cff +X) '^i' 
model is preferred by the CMB-|-BAO-|-i?o combination, 
with Ax^ = 7.9. As we will discuss later in this paper in 
§ |9.2| and Figure [17) this model space also improves the 
consistency between the CMB, BAO, and Hq datasets. 

The ACDM-l-w model is not significantly preferred by 
either the CMB or the CMB+BAO+i7o combination. 
Additionally, we find that adding the SPT bandpowers 
do not significantly improve w constraints. Therefore, we 
do not present constraints on ACDM-t-w in this paper, 
although we sometimes include w when exploring the 
effect of parameter degeneracies on the constraints for 
other parameters. 

Next we explore the second question: how do these ex- 
tensions affect the aforementioned tension between the 
CMB and low-redshift (BAO and Hq) measurements? 
This question is most clearly explored in the context of 
Figure [3] All the panels of this figure have the same x- 
and y-axes as Figure [2) The B AOboss dataset has been 
singled out for two reasons. First, it is the most precise 
low-redshift BAO measurement; second, it has the most 
tension with the CMB and Hq data in the ACDM model. 
Each panel of the figure shows color-coded samples from 
a CMB-only MCMC for a different cosmological model. 
The color coding reflects values of the extension parame- 
ter, thus showing how the chosen parameter moves within 
the CMB constraints in this parameter space. 

Among the 1-parameter extensions, ACDM-fiVoff is 
most effective at reconciling the CMB, BAO, and Hq 
datasets. A thorough discussion of the ph ysica l and 
observational effects of A^eff is reserved for § 18. how- 
ever, we briefly preview the important aspects here. The 
expansion of the CMB likelihood volume towards the 
B AOboss +^^0 constraint arises because increasing Ncs 
increases the expansion rate at early times which reduces 
the sound horizon length r^. To match the measured 
CMB acoustic peak locations, one must also decreas e Da 
by increasin g Hq. The net result, as discussed in |Hou 



et al. (2011), is an increase in Hq with minimal change 
to Ts/Dv- 'i'hus, the contours for the ACDM-fiVcfr case 
resemble those for ACDM except expanded horizontally 
in Hq. 

The extensions which become important primarily at 
late-times (Hfe, w, and '^m^) have a smaller effect on 
Ts and instead influence the late-time geometry and thus 
the inverse distance measures l/Dy and Hq. As a result, 
these extensions stretch the CMB contours along lines 
with positive slope in the HQ-Vg/Dy plane, as seen in 
the right column of Figure |3] Low-redshift datasets are 
especially important for constraining these extensions. 

In summary, the CMB data prefer the dris/dhik and 
Yp extensions to ACDM because these extensions al- 
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Fig. 4. — The impact of the angular diameter distance to last 
scattering, Da, and the lensing information on the curvature con- 
straint from the CMB. The contours show the 68% and 95% C.L. 
contours when lensing information is included. The colored points 
are samples from the posterior distribution after marginalzing over 
, which effectively removes all lensing information. The allowed 
range of Cf^Q values remains finite in the absence of lensing in- 
formation, because the CMB still places (weaker) constraints on 
the ACDM-f f2fe model parameters. The color reflects the curva- 
ture value as indicated by the color bar on the right side of the 
plot. The lensing sensitivity of the data is clearly crucial to the 
curvature constraint from the CMB data. 

low for an increasingly red tilt at higher multipoles. 
The preference for these two extensions remains in the 
CMB-l-BAO-l-iJo combination. Although the quality 
of fit improves less for the ACDM-|-iVoff model, this 
model maximizes the consistency between the three 
datasets. The CMB-t-BAO+i/o data also prefer the 
ACDM+A^off +X] '^i' model, although this preference is 
less significant (given the 1 extra dof) than that for 
running. The three datasets are completely consistent 
within this two-parameter extension. 

5. CURVATURE 
The SPT-f WMAP7 constraint on the mean curva- 



tureof the universe has been presented by [Story et al.| 
(|2012|. In brief, the addition of the SPT bandpowers 



tightens the WMAP7 constraint by a factor of five to 
flk = "O-OOS^Q Q^^g. Here we briefly present the physical 
underpinnings of how the SPT bandpowers tighten the 
constraint on fl^. 

The constraint on curvature tightens due to the sensi- 
tivity of the SPT bandpowers to low-redshift information 
through gravitational lensing. S12 used the combina- 
tion of SPT and WMAP7 bandpowers to report an 8.1 cr 
detection of gravitational lensing of the CMB. Without 
the lensing information, the CMB constraints exhibit a 
strong degeneracy between curvature and the dark en- 
ergy density. In essence, the primary CMB anisotropy 
is exquisitely sensitive to the angular-diameter distance 
to last scattering, but this is only one number. There- 
fore, if there are n parameters that only affect late-time 
geometry and structure, the CMB data will only tell us 
about one direction in this n-dimensional space. This is 
insufficient to distinguish between the two late-time pa- 
rameters flk and Q\. Lensing adds a second late-time 
measurement, breaking this degeneracy. 
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TABLE 2 

Improvements to the quality of fit for extensions to the ACDM model 





ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


ACDM 


Datasct 


+dns/dlnk 




+Yp 






+UI 


+ ^^011 + 1] 




CMB 


4.9 


1.1 


4.4 


2.4 


0.3 


0.0 


2.6 


5.5 


CMB+BAO+Ho 


7.4 


3.4 


5.3 


2.8 


2.0 


0.2 


7.9 


5.4 



This table shows how extensions to the ACDM model change the quality of the fit to the data for the CMB 
and the datasets combination of CMB+BAO+-ffo- For each dataset we report Ax^ defined as the reduction 
in from the best-fit ACDM model to the best-fit ACDM+(extension) model. A large Ax^ relative to the 
baseline ACDM model indicates that the extended model is a much better fit to the data. For instance, if 
Ax^ = 4 (6.2) for a given extension with one (two) degrees of freedom, this means that the data favor this 
extension at 2 a. Of the one-parameter model extensions considered here, both datasets have the strongest 
preference for non-zero running. Of the two-parameter model extensions considered here, the combined data 
prefer the ACDM+A^cff+5Z "^i' model most strongly. 



We demonstrate the effect of including lensing on cur- 
vature constraints from the CMB in Figure |4j The axes 
in the figure are the two late-time quantities constrained 
by the CMB data: the angular-diameter distance to last 
scattering, Da^ and the angular power spectrum of the 
lensing potential evaluated at I = 100, Cf^Q. The color- 
coded points are drawn from an SPT -t- WMAP7 MCMC 
chain with lensing information removed by marginaliza- 
tion over A^.^^ We can see the impact of the geometric 
constraint, confining the samples to a narrow range of 
Da values. Also, despite removing all lensing informa- 
tion, we see that the ACDM+rifc model predicts a finite 
range of values of the lensing power. 

As Figure |4] shows, varying flk leads to significant vari- 
ation in the CfgQ direction in the Da — Cf^Q plane. The 
direction of the response in this plane is an important as- 
pect of why lensing plays a significant role in constraining 
ftk whereas it does not for 'Y^m^. The effect of lensing 
information (i.e., requiring = 1) is shown by the black 
contours. The observed lensing amplitude rules out the 
negative curvature tail, and as mentioned above, tight- 
ens the curvature constraint from the CMB by a factor 
of five. 

Better constraints on curvature are possible by includ- 
ing low-redshift probes. For instance, CMB-t-BAO+iJo 
leads to a constraint on curvature of $7^ = —0.0061 ± 
0.0040. Even with low-redshift datasets included, SPT 
data remains important. Without the SPT bandpow- 
ers (and lensing information therein), the uncertainty 
would be roughly 15% larger: = -0.0019 ± 0.0047 
for WMAPI+^KO+Hq. 

6. MASSIVE NEUTRINOS 

We now consider allowing non-zero sum of neutrino 
masses. As we saw in § |4j the CMB data provide only 
a weak preference for this extension. However, once we 
include datasets sensitive to late times, the support for 
non-zero neutrino mass is fairly robust to the particular 
choice of the datasets and additional extensions of the 
ACDM model. 

The range of possible neutrino masses is currently con- 
strained from both above and below. Neutrino oscil- 
lation experiments place a lower bo und on the sum of 
the neutrino masses of > 0.058 eV (Hamish Robertson 



^ The current CMB data are insensitive to the small shape 
variations in cf"^ and hence marginalizing over a scaling parameter 
removes all significant lensing information. 



20091. Neutrino oscillations are insensitive to the sum 
ot the neutrino masses above this minimum. The best 
laboratory upper bound on any single neutrino mass is 
0.3 — 0.6 eV (90% CL) depending on the nuclear ma- 
trix elements adopted for the electron neutrin o from the 
Kamland-Zen double beta decay experiment ( Gando, et 
2012)). Therefore, the sum of the neutrino masses 



al 



must lie between the limits 0.06 eV < < 1.8 eV. 

Cosmological observations can provide significantly 
stronger upper limits on the sum of the neutrino masses. 
The strongest constraints on ^ rrii, come from combin- 
ing CMB measurements with low-redshift information. 
In this section, we include low-redshift measurements 
of the BAO feature. Ho, the halo power spectrum de- 
rived from SPSS luminous red galaxies (LRGs) from 



Reid et al. (2010), and the cluster abundances from the 
SP'l' SZ-selected cluster catalog (SPTcl) described by 
Reichardt et~al] ( |2012b[ ). The LRG, Hq, and WMAP7 
data have been used previously to set an upper limit 
on the neutrino masse s of < 0.44 eV (95% CL) 

(IKomatsu et al. 2011). A somewhat tighter limit of 

17'' 



201 ID . 

< 0.32 eV (95% CL) has been placed us ing SZ- 

selected gala xy clusters, CMB, BAO, and Hq by IBenson] 



et al. (2011). Similar limits of about 0.3 eV come from 
comb ining CMB-|-BAO-|-J?n with galaxy clustering d at a 



(Thomas et al. 2010 Riemer-S0rensen et al. 2012a de 
Putter et al.|2012| or A-ray cluster abundance and clus 
ter gas mass tra cnon measurements ( Mantz et al.pOlO ) 

n 



As we note in § |6.2[ the inclusion of the new BAD data 
(particularly the BOSS point) increases the preference 
for nonzero mass. 

In this section, we present updated constraints on the 
sum of the neutrino masses incorporating the new SPT 
bandpowers. We first discuss how measurements of the 
CMB, structure growth, and geometry constrain massive 
neutrinos. We then present constraints for combinations 
of the SPT bandpowers and other datasets. We finally 
consider potential degeneracies between the sum of the 
neutrino masses and other parameter extensions such as 
running of the spectral index. 

6.1. Cosmological effects of massive neutrinos 

In this subsection, we give a physical description of the 
effects of the massive neutrinos on cosmology. We start 
with the effects massive neutrinos have on the CMB. We 
note that with the current experimental precision, infor- 
mation from CMB lensing is not important to the CMB 
constraints on ^ rui,. We then move on to consider other 
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Fig. 5. — This figure shows the effects of massive neutrinos on 
the CMB power spectrum. The curves show the ratio of mod- 
els with X] mi/ = 0.5 (dashed green) and 1.0 eV (solid blue) 
respectively to the best-fit ACDM (J^m^=OeV) model spectrum 
for SPT-j- WMAP7. When increasing m^y we adjust fi^ down- 
ward to keep 8s fixed. On large scales, we see a reduction in the 
power added by the late-time ISW effect. On intermediate scales 
below the neutrino free-streaming length, we see a reduction in the 
power contributed by the early ISW effect. On scales smaller than 
the neutrino free-streaming length, the more rapid decay of grav- 
itational potentials boosts the early ISW power. The amplitude 
of the early ISW effect is damped at Z > 500 by averaging over 
multiple positive and negative contributions. 

effects in the low-redshift universe. 

To understand neutrino mass constraints from CMB 
data, we must understand how the predicted CMB power 
spectrum changes with neutrino mass. In the standard 
thermal history of the universe, massless neutrinos have 
a temperature corresponding to ~ O.lTeV at the epoch 
of last scattering. The scale at which masses start to 
have an appreciable effect is set by this temperature to 
be ^rrii, « 3 x 0.17 eV.^^ Neutrino masses well below 
this value have no imp act on primary CMB anis otropy. 
Hu fc Dodeison| ( |2002[ ) and |Ichikawa et al.| ( [20051 ) study 
in detail the imipact of higher masses on the CMB and 
find the dominant impact is due to the ISW effect. 

In a matter-dominated universe with zero mean cur- 
vature, gravitational potentials remain constant to first 
order in linear perturbation theory. Adding components 
that do not cluster, while keeping the curvature fixed 
to zero, increases the expansion rate which causes the 
gravitational potentials to decay. As photons traverse 
these decaying potentials on their way toward the ob- 
server, new anisotropics are created by what is called 
the Integrated Sachs-Wolfe (ISW) effect. The ISW 
anisotropy is generated in the ACDM model both at 
early times, as photons free stream immediately after 
decoupling through a not-completely-matter-dominated 
universe (the early ISW effect) and at late times after 
the cosmological constant becomes important (the late 
ISW effect). 

We illustrate how the ISW effect changes with neutrino 
mass in Figure [5] In this figure, we plot the ratio of 
Ci at either ^TOy= 0.5 or 1.0 eV relative to a fiducial 
ACDM-hX; model Cf"^ with ^ = 0.0. The baryon 
density ujh, cold dark matter density lOc, and the sound 
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For simplicity we assume three families of neutrinos with de- 
generate masses. 



Fig. 6. — This figure illustrates the degeneracy between rig and 
J^m^, and its role in the SPT+ WMAP7 preference for nonzero 
neutrino masses. The contours are the 68% and 95% confidence 
intervals in the ACDM-fJ^ m,, parameter space for WMAP7 (red) 
and SPT+ WMAP7 (blue). The SPT data prefer a lower value 
of ris than WMAP7, which leads the CMB data to prefer higher 
T,m.,y. 

horizon scale 9s are fixed between the three models - 
^ TO^ and Q\ vary. Three regimes are labeled in the 
figure: a reduction of power due to the late-time ISW 
effect at Z < 20, a reduction of power due to the early 
ISW effect at 20 < Z < 100, and an increase in power 
due to the early ISW effect at 100 < Z < 500. We briefly 
explain these three regimes in the next paragraphs. 

As ^ nil, increases with w;, + Wc fixed, the expansion 
rate increases at early times. Therefore, J^a must de- 
crease (increasing Da) to keep 9s fixed. Without this 
adjustment to JIa, 9s would change, primarily due to the 
change in Da- With this adjustment, we find that in 
the mass range of interest, H{z) increases relative to the 
^TO^ = model at z > 1 and decreases at z < 1. The 
decreased expansion rate at 2; < 1 results in less decay 
of the gravitational potential on very large scales, and 
therefore a reduction in the contribution to the power 
from the late-time ISW effect. The net effect is less power 
at I < 20. However, the large cosmic variance at these 
low multipoles makes the CMB data largely insensitive 
to the reduced power. 

On scales shorter than the neutrino free-streaming 
length, the increased expansion rate just after photon 
decoupling enhances the decay of gravitational poten- 
tials and thus enhances the early ISW effect. On scales 
longer than the free-streaming length, the early ISW ef- 
fect is suppressed; the clustering of neutrinos prevents 
the potential from decaying more rapidly and, due to the 
increased expansion rate, there is less time for the early 
ISW effect to accumulate. The dividing line in multipole 
space between these two regimes increases with '^m,^. 
The magnitude of the ISW effect decreases with increas- 
ing I as cancellations between an increasing number of 
positive and negative contributions washes out the sig- 
nal, becoming negligible by I ~ 500. 

The reduction of power at I < 100 and increase of 
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power at I > 100 shown in Figure [5] can be mimicked 
by an increase in ojf, and n^. As a result, ^ rrii, is anti- 
correlated with lui, and Ug- Figure [6] shows this anti- 
correlation between ^ and ris . 

Although the ISW effect is limited to angular scales 
covered by WMAP, SPT bandpowers contribute signif- 
icantly to the neutrino mass constraints by alleviating 
these degeneracies between ^TO;^, Us, and Uf,. In par- 
ticular, SPT alone prefers a lower value of Us relative to 
WMAP7 (see Figure[l]), which causes the preferred value 
of ^ m^, to increase when SPT data are combined with 
WMAP7. This shift in the SPT-preferred value of dis- 
appears with a freely varying Yp or dus/dhik, thus the 
shift in for the combination of SPT-K WMAPl dis- 

appears as well when these parameters are added. These 
shifts in the preferred value of ^ nii, are seen clearly in 
Figure [7j 

Next, we examine the other observable consequences 
that massive neutrinos have on the geometry and the 
growth of structure at late times. The decreased ex- 
pansion rate aX z < 1 directly affects the late-time 
geometry, and geometric observables. BAO measure- 
ments constrain the late-time expansion rate through 
Tsizti^eLg) I Dv{,z). Although increasing '^my from zero 
to 1.5 eV has some effect on ^^(zdrag), the dominant ef- 
fect is the increase in Dv{z) as VLa decreases. There- 
fore rs{zdj:a.g) / Dv{z) and Hq respond similarly, decreas- 
ing with increasing neutrino mass as shown in Figure |8] 

For z > 1, the increased expansion rate suppresses 
growth on scales below the neutrino free-streaming 
length (^140Mpc for a 100 meV neutrino). On larger 
scales the neutrinos can cluster, which counteracts the 
suppression. For z < 1, the reduced expansion rate 
boosts the rate of structure growth on all scales. We has- 
ten to point out, though, that the net effect on integrated 
growth to z = on small scales is that of suppression. 

In summary, the summed neutrino mass is constrained 
by the CMB primarily through the early ISW effect. The 
SPT data improve ^ constraints indirectly via an 
improved determination of Ug and the anti-correlation of 
^ mil, with Us . Adding probes of low-redshift geome- 
try and structure growth further improve constraints on 
^ nrii, . We discuss these constraints in the next section. 

6.2. Cosmological constraints on massive neutrinos 

We now look at the cosmological constraints on neu- 
trino mass. In order to discuss changes in the poste- 
rior probability density for ^ near the prior at zero, 
we will sometimes quote a different confidence interval 
than the rest of the paper. In Table [3] and this sec- 
tion, we report the estimated confidence interval ^ to^ € 
[xi,x2] (68% CL) such that 



P(^m,)d(^m,) = 0.68, 



(5) 



where P{xi) = P{x2) and P(Ym^) is the normalized 
neutrino mass posterior probability density. When we 
find P(0) > P{x2), we report the confidence interval as 
[0, X2] ; this is an upper limit. We will also report the peak 
of P(^my). For cases with a non-zero 95% confidence 
lower limit as defined above, we will continue to report 
the median and 1 a err or w hich is used in the rest of the 
paper and defined in § 2.3 
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Fig. 7. — The constraints on from the CMB data de- 

grade significantly when we decouple the tilt of the primordial 
power spectrum on small and large scales. The marginalized one- 
dimensional posteriors of Y^rrii, from WMAP7 {SPT+ WMAP7) 
data are shown with the orange solid (black solid) line. The SPT 
data prefer a lower value of Us than WMAP7, which leads the 
CMB data to prefer higher "^m^. If the tilt of the primordial 
power law on small and large scales is decoupled, whether directly 
by introducing running (blue dashed line) or indirectly by freeing 
the damping scale with the introduction of Yp (blue dot-dashed 
line), the shift towards higher neutrino masses in the CMB likeli- 
hood is substantially reduced. 

The SPT data prefer lower values of Ug than those 
preferred by WMAP7, which translates into a preference 
for higher neutrino masses. For the SPT-I- WMAP7 data, 
we find the position of the peak moves to 0.93 eV with 
a 68% confidence interval of "i^ S [0.41, 1.34] eV. The 
95% CL upper limit from the CMB iaJ2m^ < 1-60 eV. 

We now consider low-redshift probes of geometry: 
BAO and Hq. As mentioned above, these low-redshift 
geometric probes can strengthen inferences of neutrino 
mass through the impact on the expansion rate. With 
massive neutrinos, the CMB data is compatible with ei- 
ther the BAO or Hq measurements (see Figures [3] and 
^ , though the resulting constraints on ^ m^^ are quite 
oifferent. Adding BAO data to the CMB data tightens 
the neutrino mass constraint significantly; for the com- 
bination of CMB-I-BAO, we obtain a 2.4 ct preference for 
nonzero neutrino masses with ^ — (0.49 ± 0.20) eV. 
The 95% confidence interval is S [0.11, 0.88] eV. 

On the other hand, the Ho dataset prefers higher val- 
ues of Hq, corresponding to lower values of '^m^. The 
likelihood for CMB-f-i/p peaks near zero neutrino mass, 
with an upper limit ofJ2^v< 0.48 eV (95% CL). Com- 
bining all three datasets produces a neutrino mass con- 
straint that lies between the CMB-I-BAO and CMB-|-iJo 
constraints. For the combination of CMB-fBAO-l-i/o: 
we find the posterior peaks at 0.33 eV with a 68% con- 
fidence interval of "i^ e [0.16, 0.51] eV. The 95% CL 
upper limit is J2''^f < 0.66 eV. The constraints from 
these combinations are shown in Figure [9j As an aside, 
we note that as in the CMB-only case, the small-scale in- 
formation from SPT data is a significant contributor to 
the suggestion of nonzero neutrino masses. Without the 
SPT bandpowers, the neutrino mass probability density 
for WMAP7+JiAO+Ho peaks just above zero mass with 
a 95% CL upper hmit of ^ < 0.48 eV. 

We next turn to low-redshift probes of structure 
growth: LRGs and galaxy clusters. For the CMB-I-LRG 
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TABLE 3 

ACDM+J]m,^ RESULTS FROM VARIOUS COMBINATIONS OF DATASETS 



Datascts 






m,u [cV] 




1 n2 









Us 




^0 


''"8 


Combinations 


68% 


CL 


95% CL 


Peak 














[km Mpc^i 




CMB 


[0.41, 


1.34] 


[0,1.60] 


0.93 


2.184 


± 


0.041 


0.948 


± 


0.014 






CMB+BAO 


[0.29, 


0.68] 


[0.11,0.88] 


0.48 


2.207 


± 


0.034 


0.957 


± 


0.008 


67.1 ± 1.2 


0.712 ±0.048 


CMB+ffo 


[0,0 


.24] 


[0, 0.48] 


0.01 


2.233 


± 


0.035 


0.965 


± 


0.009 


71.6 ± 1.8 


n 7KO+0.031 


CMB+LRG 


[0.17, 


0.61] 


[0, 0.78] 


0.39 


2.204 


± 


0.036 


0.956 


± 


0.009 


66.4 ± 2.5 


0.731 ± 0.048 


CMB+SPTcL 


[0.28, 


0.80] 


[0, 1.03] 


0.55 


2.196 


± 


0.040 


0.954 


± 


0.011 


65.ll|i 




CMB+BAO+i^o 


[0.16, 


0.51] 


[0, 0.66] 


0.33 


2.219 


± 


0.034 


0.959 


± 


0.008 


68.4 ± 1.0 


0.741 ± 0.046 


CMB+BAO+//0+ SPTcL 


[0.21, 


0.43] 


[0.10,0.54] 


0.32 


2.221 


± 


0.034 


0.960 


± 


0.008 


68.3 ± 1.0 


0.739 ±0.027 



This table shows the results for the ACDM+ model. The confidence level intervals shown indicate the region of highest 

probability density that contains either 68% or 95% of the probability. The 'Peak' value is the peak of the posterior distribution. 
We include some other parameters of particular interest for this extension. Note in particular that high peak values of rriv 
correspond to low values of ris , low values of Hq and low values of erg . 



dataset, we find the posterior probability peaks at 
0.39 eV with a 68% confidence interval of X^^^f ^ 
[0.17, 0.61] eV. The 95% confidence interval is an up- 
per limit at '^rriy < 0.78 eV (95% CL). Instead, adding 
the SPT galaxy cluster samplers (SPTcl) to the CMB 
leads the posterior probability to peak at slightly higher 
masses 0.55 eV with a 68% confidence interval of ^ e 
[0.28, 0.80] eV. The 95% CL upper limit is J^^i' < 
1.03 eV. Both the CMB+LRG and CMB+SPTcl com- 
binations show a preference for non-zero neutrino mass 
at just under 2 a. 

Finally, we combine the CMB, geometrical, and large- 
scale structure observations. Because the SDSS LRG and 
SDSS BAO galaxy samples overlap, we do not include the 
LRG dataset in this combination. In this case, we find 
(for CMB+BAO+Fo+SPTcl): 

= (0.32 ±0.11) eV, (6) 

^m^e [0.01, 0.63] eV (99.7% CL) . (7) 

This data combination gives a 3 cr detection of posi- 
tive neutrino masses. Without BAO, the preferred mass 
drops significantly while the uncertainties increase; the 
constraint is consistent with no mass at 68% confi- 
dence, and we end up with an upper limit of ^ m,y < 
0.39eV(95%CL). Without Hq, the likelihood shifts 
to higher masses without degrading the uncertainties, 
(0.40 ± 0.11) eV. 
The are several drivers behind the ^ constraint 
(Equation ro| derived from the dataset combination of 
CMB-hBAO+Ho+SPTcL: 

• CMB: As discussed above, the CMB data con- 
tribute to this constraint by enforcing the height of 
the first and second acoustic peaks (from WMAP7 
data) and the slope of the damping tail (from SPT 
data) . The latter leads to a lower Us , with ^ rrii, 
thus increasing to compensate for the changes to 
the power in the WMAP7 region. The resulting 

Note that the interpretation of the measured galaxy number 
counts depends on the Tinker mass function (Tinker et al. 2008]|. 
Although the Tinker mass function was not originally c alculated 
for mas sive neutrinos, later papers (MaruUi et al. 2 0lT| [ichiki &[ 
|Takada||2012) have shown that if rescaled to the new value ot (Tg, 
the Tniker mass function remains accurate for massive clusters 
(M > W^'^H'^Mq) - which includes all SPT galaxy clusters. 



posterior distribution from the CMB peaks around 
J2 iT^u ^ 1 eV. These CMB constraints contribute 
~ 1 (T to the combined ^ constraint. If we 
keep only the optical depth prior from WMAP7, we 
find a confidence interval of ^ rrii, G [0.07, 0.53] eV 
(95% CL). Without the SPT bandpowers (and Us 
information) , the confidence interval is ^ rrii, € 
[0.03, 0.48] eV (95% CL). In both cases, there re- 
mains a greater than 2 a preference for nonzero 
masses. 

• BAO: The BAO measurement is an important 
driver of the combined neutrino mass constraint. 
As seen in Figure |9l the BAO measurements push 
up the preferred value of ^ nii, and tighten the con- 
straints considerably relative to the other datasets. 
The preference for massive neutrinos is common 
between the three BAO datasets. To confirm the 
robustness of the BAO preference, we have run 
three MCMC chains, dropping one of the BAO 
datasets in each chain. When dropping the SDSS, 
WiggleZ, or BOSS dataset, we find, at 95% CL, 
X;™,. G [0.10,0.55], [0.10,0.53], or [0.03, 0.47] eV 
respectively. In short, the BAO data prefer massive 
neutrinos, and the results are robust against drop- 
ping any single BAO observation. The largest shift 
of 0.6 cr is introduced by dropping the BOSS mea- 
surement. However as previously noted, if all BAO 
measurements are removed and Hq information is 
included, the preference for massive neutrinos dis- 
appears. 

• SPTcl: Adding the SPTcl galaxy cluster sample 
reduces the neutrino mass uncertainties by a fac- 
tor of 1.6 without significantly changing the me- 
dian value, as can be seen by comparing the value 
without clusters - = (0.34 ± 0.18) eV - to 
the value with clusters - ^ m„ — (0.32 ± 0.11) eV. 
In principle, this result depends on the accuracy 
of the cluster mass calibration, which was deter- 
mined through an X-ray based scaling relation. To 
test this dependency, we have doubled the uncer- 
tainty on the normalization of the X-ray scaling 
relation and find little change in the distribution, 
X; TO,. = (0.34±0.12) eV. We conclude that the re- 
duction in uncertainties is insensitive to the specific 
prior on the mass scaling relation. 
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Fig. 8. — This figure shows how low-redshift datasets aff'ect 
the neutrino mass constraints. In the planes of ^'^'^ three 

quantities inferred at low redshifts, we show the marginalized 
two-dimensional posteriors from different combinations of datasets 
shown by the legend in the central panel. Top: constraints are 
plotted in the "^rrii, — rs(zdrag)/-Dv(0.35) plane with the SDSS 
BAO data point shown by the grey band as the 1 cr region. Middle: 
similar to the left panel but showing rs{z^^i^g)/ Dy{0.57) with the 
BOSS BAO data point for the gray band as the 1 a region. Bot- 
tom: similar to the top and middle panels but showing Hq with the 
grey band showing the 1 a region. The tension between BAO and 
Ho can be seen in the three panels. The CMB+BAO preference 
for massive neutrinos is clearly visible. 



We conclude that the preference for nonzero neutrino 
masses is coming from three independent sources and 
remains at > 2 cr significance as fong as we keep cluster 
abundances and at least one of the two most precise BAO 
measurements. 

6.3. Degeneracies with other extensions 

A final question is to what extent the neutrino mass 
constraint is weakened by introducing additional free pa- 
rameters to the 7-parameter ACDM-f ^ model. To 
address this, we show the marginalized posteriors for the 
sum of the neutrino masses in 8-parameter models in Fig- 
ure [TO] We consider several possible additions: Yp, Ncs, 
dris/dlnk, H.^, or w. All curves are for the combination 
of CMB-fBAO-Hifo+SPTcL- 

Of these five extensions, two {dris/dhik and Qk) de- 
crease the preference for positive ^ m^, two {w and A^off) 
increase the preference, and one (Yp) has minimal effect. 
Allowing both curvature and massive neutrinos signifi- 
cantly increases the uncertainties on both parameters. 
The peak of the neutrino mass likelihood remains nearly 
unchanged at 0.34 eV with the 95% confidence interval 
expanded to ^ [0, 0.70]eV. The curvature is con- 

sistent with zero. Allowing for running of the spectral in- 
dex reduces the preference for nonzero neutrino masses, 
in this case by shifting the median down; the neutrino 
mass constraint with running is ^ nii, — (0.27±0.11) eV. 
On the other hand, allowing the dark energy equation 
of state or number of neutrino species to vary increases 
the significance of the neutrino mass detection. The 



sum med neutrino mass is c orrelated with Ncs (see § 9.2 
and Ichikawa et al. 2005) and anti-correlated with w 
(Zhao et al. 2012). treemg either of these parameters 
increases the median masses, although the uncertainties 
also increase. For a wCDM-|-^ m^, model, the mass con- 
straint is X] = (0.51 ± 0.14) eV, 3.7 cr above zero. For 
a ACDM-f A'off +X] cosmology, the mass constraint is 
J2iTiu — (0. 51 lb 0.15) eV. This last case is discussed 
in detail in §|93] We find the CMB-bBAO-b-ffo+SPTcL 
preference for neutrinos remains at greater than 95% con- 
fidence when we add additional model parameters, with 
the exception of the mean curvature of the universe. Fit- 
ting the posterior with a Gaussian, we find curvature 
reduces the preference for massive neutrinos to 1.7 cr. 

7. RUNNING OF THE SPECTRAL INDEX 

The difference in preferred Ug from SPT and WMAP7 
data suggests a scale dependence in the power spectrum 
of the primordial fluctuations. To test this, we allow the 
primordial power spectrum to deviate from a pure power 
law by introducing a logarithmic de pendence on scale k, 
a so-c alled "running spectral index" ( Kosowsky & Turner 
1995|: 



ns{k) = ns(fco) + dng/d\nk\n 



(8) 



Here, fco is a defined pivot point and dns/dln k is the run- 
ning parameter. Throughout this section, we will define 
this pivot point to be k^ = 0.025 Mpc~-^, which projects 
to I ~ 350. This pivot point is chosen to decorrelate the 
uncertainties on rig and dris/dlnk in the SPT+WMAP7 
data. 
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Fig. 9. — Left panel: This panel demonstrates how the marginahzed one-dimensional posterior distribution of "^rrii, changes when a 
single external dataset is added to the CMB dataset. The CMB constraint is shown by the black solid line. The datasets added are the 
BAG data (blue dot-dashed line), Ho measurement (blue dashed line), LRG sample (black dotted line), and SPT cluster data (orange solid 
line). Right panel: This panel demonstrates how the marginalized one-dimensional posterior distribution of X]"^"^ changes when a single 
external dataset is removed from the combination of CMB-|-BAO+i?o+SPTcL- The datasets dropped are the SPT bandpowers (orange 
solid line), SPTcl (blue dashed line), BAO data (blue dot-dashed line), and Hq measurement (black dotted line). The marginalized 
posterior for the combined CMB-f BAO-|-//o+SPTcl dataset is shown by the black solid line. 
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Fig. 10. — This figure illustrates the robustness of the neutrino 
mass detection to other parameter extensions. The marginalized 
one-dimensional posteriors for are shown for two-parameter 

extensions to ACDM for the combined CMB-(-BAO-f Hq-I-SPTcl 
data sets (for ui, SNe are used instead of -ffo)- Allowing significant 
curvature or running can significantly reduce the preference for 
nonzero neutrino masses (to 1.7 and 2Acr respectively). Other 
extensions increase the preference for positive neutrino masses. 

The running parameter dris/dhik is predicted to be 
undetectable by most inflationary theories, and a de- 
tection of non-zero dns/dlnk could prov ide information 



adding the SPT bandpowers to WMAP7 dramatically 
reduces the allowed likelihood volume, and leads to a 
preference for negative running. For SPT+ WMAP7, we 
find: 

dn^/d In fc = -0.024 ±0.011. (9) 

The probability of negative running with the CMB data 
is P{dns/d\nk < 0) = 98.6%, equivalent to a 2.2a 
Gaussian preference. The preference for negative run- 
ning remains above 2 a when the beam uncertainties are 
doubled , wh en the priors on the three foreground ampli- 
tudes (§ 2.1 1 are removed, and when we restrict the SPT 
bandpowers to the I < 1500 range; we conclude that this 
preference does not depend strongly on the SPT experi- 
mental beam uncertainty or foreground modeling. 

Adding BAO and Hq data marginally improves the 
constraints and shifts the median to more negative 
values. The combination CMB-fBAO+i?o constrains 
drig/dhik to be 



dn^/d \nk = -0.028 ± 0.010. 



(10) 



running. 
dns/d lnk 



Komatsu et al 
0.020 (95% CL) 



< 



about the inflationary potential (Kosowsky fc Turner 
19951, or point to models other than mhation. i'here 
have been a number of rec ent C MB constraints on 

(2011|) obtain - 0.084 < 
WMAP7. [Dunkley 
et al.| ( |2011[ ) find dnjdlnk = -0.034 ± 0.018 from' 
WMAP7+ACT. Kll use the combination of WMAP7 
data and the first 790 deg^ SPT survey data to obtain 
dus/dlnk = — 0.024 ± 0.013, a preference for negative 
spectral running at 1.8 ct. 

7.1. Constraints on dug/dhik 

We now look at the constraints on dus/dXn k from the 
SPT bandpowers. As shown in the left panel of Figure 11 



The small shift from the CMB constraint is driven al- 
most entirely by the BAO data, which prefers smaller 
values of H^. The probability of negative running with 
the CMB+BAO+iJo data is P{dnjd\nk < 0) = 99.7%, 
equivalent to a 2.7 a Gaussian preference. 

We next consider how robust the preference for neg- 
ative running is to other model extensions. Introduc- 
ing tensor perturbations (by making the tensor-to-scalar 
ratio r a free parameter) is a particularly natural ex- 
tension in the context of inflationary models. How- 
ever, the preference for running cannot be reduced by 
non-zero tensor perturbations. Tensor modes would in- 
crease the power at large scales without affecting small 
and intermediate scales. Adding tensor modes, there- 
fore, increases the preference for negative running. In 
a AGDM+r+dus/dlnk cosmology, we obtain dUs/dlnk 
= -0.045 ± 0.016 for the CMB data and dnjdlnk 
= -0.046 ± 0.015 for the CMB-hBAO-hiJo- 

More generally, the preference for negative running 
might be reduced by any extension that also effec- 
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Fig. 11. — The data prefer negative running; this preference remains when additional parameters are freed. We show the marginal- 
ized two-dimensional posteriors for ris (at k = 0.025 Mpc~^) and dus/dlnk. Left panel: the Us- dus / d\nk joint distribution in the 
ACDM-fdns/dln fc model for different combinations of the datasets. Right panel: the n^- dUs/dlnk joint distribution obtained by 
marginalizing A^j.ff (red) and Yp (blue) for WMj4P7+SPT data. The purple filled contours show t he cas e with Yp marginalized with 
the prior Yp < 0.294, the 2 a upper limit of the solar initial helium abundance by Sercnelli & Basu] |2010| . We show the region where 
|dns/dlnA;| < (1 — Us)^ in gray. Inflation models with slow-roll expansion that can be terminated at second order make predictions in this 
region. 

certainties, dnJdXnk = -0.034 ± 0.017. We find the 
largest increase in the uncertainties for dng/dlnk +5^,, 
dus/dlnk = -0.020 ± 0.024. Therefore, the significance 
of the preference for running is most reduced by freeing 
Yp and, to a lesser extent, N^^g, and we plot constraints 
for these extensions in the right panel of Fig. [TT] as well 
as Fig. [T2| 

As wilTbe discussed in § [Sl the primary effect of Yp and 
Ncff is on the damping scale - which obviously mimics 
running. This degeneracy expands the uncertainty on 
dus/dhik. The shift in the ISW effect due to massive 
neutrinos does not mimic the scale dependence of run- 
ning very well and has minimal effect. Freeing A'efi and 
l^yield different results because, as will be discussed in 
§ |9l A'eff also changes the locations of the acoustic peaks 
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Fig. 12.— The preference oftheSPT-|-VFM^P7 data for negative 
running remains as additional parameters are freed. The marginal- 
ized one-dimensional posteriors for dns/dln k from the CMB in the 
ACDM-|-dns/dln fc model is shown by the black solid line. The 
preference for negative running is reduced when Yp is allowed to 
vary (blue dashed line), although this depends on extremely high 
helium abundances. If an upper limit based on solar abundances 
is set on Yp, the preference for negative running is mostly restored 
(blue dot-dashed line). Allowing Afefi to vary (orange solid line) 
increases the preference for negative running. 



tively allows for an increasingly red, scale-dependent 
tilt. We thus consider how the running constraint 
changes for the CMB data alone when freeing ^ m,^ , 
A'eff, or Yp. For dus/dhik -|-^m^, we find dns/dlnk = 
—0.020 ± 0.012, similar to the ^ rrii, = eV constraint 
of dnjdlnk = -0.024 ± 0.011. For dn,/dlnk +Ncs, 
we find a shift to more negative running and larger un- 



- which running and Yp do not. This reduces the de- 
generacy between the parameters and also explains the 
shift in the preferred value for running. With dng/dlnk 
+Ncg free, the model is free to move to lower values of 
Ncff to better match the observed peak locations while 
compensating for the decreased damping with negative 
running. 

Among the models considered, the preference for nega- 
tive running is only removed by freeing the helium abun- 
dance. However, the inferred helium abundance is in 
tension with the 2 a upper limit of the protosolar mea- 
surement that will be discussed in ^ |8.2 Forcing Yp to be 



less than this 2 a upper limit of Yp < 0.294 degrades the 
effectiveness of Yp in removing the preference for negative 
running. This point is illustrated by the purple contours 
in the right panel of Figure [TT] and the blue dot-dashed 



line in Figure 12 



Excluding values of Yp above 0.294, 
the preference "Tot running is robust against other model 
extensions. 



7.2. Implications for Inflationary Models 
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We now explore the implications of the constraint on 
running for models of inflation. Although we expect Us 
to have some scale dependence, dns/dlnk is predicted 
to be undetectably small with current data in single- 
field, slow-roll inflation models. These models predict 
dus/dlnk to be second order in the slow roll parame- 
ters, and thus of order (1 — n^)^ unless the poten tial ex- 
periences signi ficant jerk in th e observable range ( Chung 
^ al. 2003; Fin elh et al.|2010| ). A significant detection of 
running is a potential problem for both small and large- 
field inflation models. Small-field infiation models pre- 
dict r ~ 0, thus one can interpret the ns-dng/dlnk con- 
straints shown in Figure 111] as a direct test of small-field 
models. The observed value of dug/d In k disfavors single- 
field small-field inflation with negligible jerk at the ^ 2a 
level (- 2.7 cr for CMB+BAO+ffo)- Large-field inflation 
models, which predict non-negligible r (and thus more 
negative running), are even more disfavored. 

Models with significant jerk at observable scales can 
allow for larger absolute values of running f 
Peiris;2006), however for rig < 1, these mode 
diet a number of e-foldings N < 30. '^'^ Here e" gives 
the increase in the scale factor between the time when 
the observable scale leaves the horizon and the end of 
inflation. Such small values of N are incompatible with 
a standard cosmological history (Liddle & Leach 2003), 
as they imply an energy scale of inflati on which is be- 
low t he electroweak phase transition (Easther & Peiris 
2006 1 . We therefore treat the prediction that dUs/dlnk 
is of order (1 — ris)^ as a fairly robust one for single-field 
slow-roll models. 

Finally, we mention that the choice of a constant 
drig/dlnk as a departure from power law behavior in 
the primordial spectrum is not the only possible choice. 
It is also possible that the runnin g is significan t over only 



Easther &; 
s also pre 



Cline & Hoi 



(2006) ar 



a small range of wavenumber fc; 
gue that such a parametrization fits more easily into the 
slow-roll framework and allows inflation to continue for 
a longer number of e-foldings. 

As we conclude this section on running, it is worth re- 
iterating that the statistical preference for running is not 
greater than 3 a. However, the other model extensions 
we consider do not significantly weaken this preference 
for running. If the constraints from future experiments 
confirm this preferred value of running, it would repre- 
sents an important clue about the earliest epochs of the 
universe. 



8. CONSTRAINTS ON OTHER ONE-PARAMETER 
EXTENSIONS 

In this section, we consider other physically motivated 
models that can affect the scale-dependence of the tilt 
of the CMB power spectrum. Specifically, we test the 
effect of (1) varying the effective number of neutrino 
species (A^^sg), a probe of the standard model of particle 
physics, and (2) an increase in the primordial fraction 
of baryonic mass in helium {Yp), a possible signature of 
non-standard BBN. These extensions can increase the 
predicted damping due to photon diffusion relative to 

Allowing for fourth or higher derivatives can increase the num- 
ber of e-foldings, as is th e case with lar ge-field chaotic models with 
sinusoidal modulations ( [Kobayashi fc Takahashi 201l| or extra- 
dimensional natural inflation models ( [l^'cng et al.^2UU3^ ^ 



ACDM, which preferentially decreases small-scale power. 
The data therefore favor these extensions. Conversely, 
the data will disfavor extensions that exacerbate the dis- 
crepancy between the observed and predicted small-scale 
power. Examples of this type of extension include the 
tensor-scalar ratio r, which as discussed in S12 leads to 
larger values of n^, and the effect of early dark energy, 
which we consider below. 



8.1. 



off 



In the standard thermal history, the radiation density 
in the early Universe is given by 



Pr = + = (1 + 0.227 A^offOp-, 



(11) 



where is the photon density, a quantity measure d ex- 
tremely well by COBE/YlRkS ( |Fixsen et aLlpOe] ), and 
A"off is an effective number of neutrino species. 'I'he fac- 
tor of 0.227 is calculated assuming all the entropy of e"*", 
e~ annihilation goes into the photons. This assumption 
is not perfect; the small correction due to the neutri- 
nos gaining some entropy has traditionally been incorpo- 
rated into NcE- For the sta ndard model of thr e e neu- 
trino species, A^eff = 3.046 (iDicus et al. 1982 Lopez 
et al.||1999] [Mangano et aLpOOS) . Although generally 
discussed in the context of extra neutrino species, A'off 
parametrizes the total non-photon contribution to the 
radiation density and thus includes contributions from 
any weakly/non- interacting, relativistic particles 



As describe d by ,Bashinsky k, Seljak^ (2004| ) and Hon 



et al.| (|2011j), the t3^-\:+WMAF data constrain Xff 
through its effects on the expansion rate at early times. 
Both the sound horizon Vg and the square of the damp- 
ing length depend on an integral with an integrand 
proportional to \/H{z). However, the angular power 
spectrum is sensitive to the corresponding angular scales 
Qd = ^d/DA and Og — rs/DA, where Da is the angular- 
diameter distance to recombination. Since the angular- 
diameter distance is not known, measurements of 9d or 
alone do not constrain the expansion rate well. The 
dependence on the unknown angular-diameter distance 
can be removed by looking at the ratio Vd/rg = Od/Og. 
As the effect A"cff has on the expansion rate is roughly 
redshift independent, its effect can be pulled out of the 
integral, so in this extension 6d/9s oc H^^^. A second 
effect of A^cff is on the phase of the acoustic oscillations. 
Constraints from this phase shift are subdominant in the 
ACDM-I-Aoff model, but become more important when 
Yp is also allowed to vary; we discuss this effect in that 
context in § |9] The Ncs constraint in the ACDM-fA^eff 
model can thus be interpreted through its effect on the 
9d/0s ratio. 

Here we hold Yp to the value predicted by BBN, which 
is a function of both N^g and cj^- These Yp predic- 
tions are implemented in CosmoMC through interpo- 
lating over a table produced using PArthENoPE vl.OO 



( Pisanti et al.| |2008|. We look at the effect of varying Yp 
independent oi BBN constraints later in this section, and 
we consider jointly varying A^off and Yp independently in 

The left panel of Figure [13| shows the CMB constraints 
on 9d/0s in the ACDM-t-Acfr model space. The parame- 
ter 0s is well-measured by both SPT and WMAP7, and 
the fractional shift in preferred 9s is small between the 
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Fig. 13. — The marginalized one-dimensional posteriors for 6d/Ss (left) and A^^ff (right) from various combinations of datasets. The 
constraint on A^cfi in the ACDM+Afcff model for different comliinations of datasets can be interpreted from the corresponding Qdl^s 
posterior. The shift in Qdl^s between broad WMAP7 (dot-dashed blue curve) and SPT (dashed blue curve) implies a preference for lower 
A^cfi with SPT than WMAP7. The N^ff posterior of SPT-I- WMAP7 is shown by the solid black curve. The tighter constraints on N^ff can 
be obtained by adding Hq (black dotted curve) and combining CMB, BAO, and Hq (solid orange curve). 



two datasets. As a result, the uncertainty in A'ofr is pri- 
marily due to the uncertainty in 0^- In ACDM, 6d is 
primarily constrained from the WMAP7 determinations 
of uJb and ujm - not by the damping tail measurement. 
Freeing Ncs greatly broadens the WMAP7 constraint on 
9d due to a degeneracy between A'cfr and oj^- With N^ff 
free, the SPT determination of 6d from the damping tail 
becomes tighter than that from WMAP7. 

Cons tra ints on A'eff are shown in the right panel of 
Figure [13] The dot- dashed blue cur v e mar ks the broad 
WMAP7 posterior; 'Komat su et al.| ( |201lD find N^s > 
2.7 (95% CL) using WMAP7 alone. As expected from 
the observed shift in 0d/Os between WMAP7 and SPT, 
the SPT data prefer lower values of N^s than WMAP7. 
Adding the SPT data to WMAP7 markedly improves 
the measurement as shown by the black curve. The joint 
SPT+ WMAP7 constraint is: 

N^e = 3.62 ± 0.48, 

representing a 20% reduction in uncertainty from the 
constraint determined from the WMAP7+K11 band- 
powers. For the CMB data, the probability that A'ofr > 
3.046 is 89%. 

We now turn to the addition of the late-time BAO and 
Hq data. When Zcq is held fixed at its well-measured 
value from the CMB, increasing the effective number of 
neutrino species results in an increased expansion rate, 
which decreases Tj. Since 9s is well constrained, this re- 
sults in a decrease in the angular-diameter distance Da, 
and thus Dy- The end result is that, for the BAO observ- 
able Ts/Dv, the changes largely cancel. The CMB+BAO 
constraint is A'eff = 3.50 ± 0.47. The direct Hq mea- 
surement is more sensitive to the value of N^ff- How- 
ever, the Hq and CMB datasets prefer similar values of 
-ffoand the preferred value of N^s hardly moves when 
Hq data is added. The resulting CMB-|-iJo constraint 
is A'eff = 3.46 ± 0.35. The significance of the prefer- 



ence for A'cff > 3.046 is largely unchanged between these 
three cases: the CMB, CMB+BAO, and CMB-^iJo pre- 
fer Ncs > 3.046 at 1.2 ct, 1.0 ti, and 1.2 cr, respectively. 

Though adding BAO or Hq data individually to the 
CMB data slightly reduces the preferred value of A'eff, 
combining all three datasets has the opposite effect of 
shifting the distribution towards slightly higher A'cfr . 
The joint CMB+BAO+Hq constraint is: 

TVeff = 3.71 ±0.35 

a 1.9(7 preference for A'cfr > 3.046. The reasons for this 
shift can be seen in Figure [3] Increasing A'cfr moves 
the CMB-predicted Hq and B^O quantities towards the 
measured values. This ability of high A'cfr to reconcile 
the CMB with BAO and Hq ha s been noted previously 



by, e.g., [Anderson et al.| ([2012|. The inclusion of the 
BAO and Hq information, leads to a small upward shift 
in the value of A'cfr; a higher value of A'cfr is required 
to reduce the tension between these datasets than that 
inferred from the CMB determination of Od/Og- 

These A'cfr constraints depend on the modeling as- 
sumptions. The effective number of neutrino species 
is partially degenerate with several other extensions 
to ACDM including dus/dlnk (which shifts the CMB 
datasets preferred value of A'cfr down to Aoff — 
2.56lg gg), the primordial helium abundance (considered 



9.1), and the summed neutrino mass (considered in 



The mild preference for A'cfr > 3 disappears when 
either dn^/dlnfe or Yp are freed; however, the preference 
remains when freeing ^ rrii, . 

In summary, CMB data primarily constrain A'cfr 
through measurements of the ratio Od/0s- Adding Hq and 
BAO measurements tightens this constraint and leads to 
a 1.9 (7 preference for more than 3 neutrino species. How- 
ever the significance of the preference is sensitive to the 
addition of further model extensions. 
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TABLE 4 

ACDM+ N^ff Results from SPT, CMB, CMB+BAO, CMB+J^o and CMB+BAO+Ho 



"As calculated from BBN theory 



Datasets 




ip 


lO^f^b/i^ 


Ho [kms-i Mpc-i] 




10 ea/Bs 


SPT+T 
CMB 
CMB+BAO 
CMB+Ho 
CMB+BAO+//0 


o qK+116 
■^■'-'"-1.03 

3.62 ± 0.48 
3.50 ± 0.47 
3.46 ± 0.35 
3.71 ± 0.35 


0.252 ± 0.014 
0.255 ± 0.006 
0.254 ± 0.006 
0.253 ±0.005 
0.256 ±0.004 


2.307 ± 0.122 
2.268 ±0.049 
2.232 ±0.045 
2.253 ±0.038 
2.247 ± 0.038 


77.ll^:? 
75.9 ±3.4 

71.4 ± 2.5 

74.5 ± 1.9 

72.6 ± 1.8 


3053 ± 192 
3136 ± 96 
3326 ± 50 
3161 ± 84 
3322 ± 49 


1.552 ±0.030 
1.562 ±0.011 
1.561 ± 0.011 
1.559 ±0.009 
1.565 ±0.009 



This table shows constraints on the ACDM+ TV^ff model. In addition to N^g, we include constraints on other 
parameters of particular interest for this extension. In this model space, Yp is related to i^i,h^ and N^g by a BBN 
consistency relation. 



8.2. Yp 

When the universe cools to T ^ 0.1 MeV, Ught nuclei 
begin to for m, a process known as big bang nucleosyn - 
thesis (BBN, [Schramm fc Turi^|1998| [Steigman||2007] . 
The primordial fraction of ba ryon ic mass in *He is de- 
noted as y^. As mentioned in § 8.1 BBN makes a precise 
prediction for the primordial helium abundance. A use - 
ful analytic form is given by Simha & Steigman ( 2008 ) , 
valid for N^ff near the standard model prediction: 

Yp = 0.2485 + 0.0016 [(273.9 loi,-6) + 100(5 - 1)] , (12) 

where 

5'2 = 1 + (7/43)(A^off - 3.046). (13) 

The factor accounts for any non-standard expansion 
rate prior to and during BBN. We use the BBN predic- 
tion for Yp in nearly all MCMCs, unless we specifically 
state that Yp is free. 

At fixed baryon density, increasing the helium fraction 
leads to increased damping. Helium has a higher binding 
energy and thus recombines earlier than hydrogen; each 
neutral helium atom effectively absorbs four free elec- 
trons at the surface of last scattering. This decreases the 
number of free electrons which increases the photon dif- 
fusion length, r^, as seen in Equation [4j Mirroring 



this single distance measure is degenerate with the angu- 
lar diameter distance to last scattering. However, since 
the sound horizon scale is nearly independent of Yp, the 
ratio of 0d/9s breaks the degeneracy with Da and allows 
the CMB data to constrain Yp independently of BBN. 

Precise experimental inferences of Yp have come from 
spectral observations of extra galactic cloud s of low- 
meta. l licity ionized gas ([ izotov et al.|2007 Peimbert ct al. 



2007 



2012 



[Izotov fc TliuaS||2010; ,Ayer et al.| [i^OlO, ,2011 
. Historically, these results have shown substantia, 
variations between groups that greatly exceed the ex- 
pected statistical variations; for a mo re complete discus- 



sion of these measurements, see e.g., Aver et al. (2010) 



Marginalizing over all possible trends ot the estimated Yp 
with metalli city increases the uncertainty substantially: 



for example, Aver et al. (2012) find Yp = 0.2534± 0.0083 



0.7 a above the BBiN prediction for Yp assuming the three 
standard neutrino species and the value of cjf, measured 
by WMAP. 

An upper limit to Yp can be inferred from protoso- 
lar abundance estimate s resulting from helioseismology; 
Serenelli & Basu (2010) find the protosolar helium abun- 
dance to be Vsoiar = 0.278 ± 0.008. This measurement 
relies upon well-studied gravitational and solar evolution 



physics to back out the protosolar abundance from cur- 
rent measurements of the helium fraction on the solar 
surface. It is difficult to imagine mechanisms by which 
the protosolar abundance is reduced below the primordial 
helium abundance. Therefore some measure of caution is 
warranted in interpreting cosmological constraints with 
primordial helium abundances above the 2 a upper limit 
ofYp< 0.294 derived from these measurements. 
After combining WMAP 7 with ACBAR and QUaD 



temperature power spectra, Komatsu et al. ( 2011 1 deter- 
mine Y„ = 0.326 ± 0.075. 'Dunkley et all POllI) obtain 



Yp = 0.313 ±0.044 by combining the ACT ana~WMAP7 

data. Using the bandpowers from the first 790 deg^ SPT 
survey and WMAP7 data, Kll find Yp = 0.296 ± 0.030. 
In this analysis, we infer Yp — 0.240 ± 0.079 with the 
SPT bandpowers alone. These results are shown in the 
right panel of Figure [14] The combined CMB data lead 
to: 



K = 0.300 ± 0.025. 



(14) 



As shown in the left panel of Figure 14, the addition 
of the WMAP 7 data breaks a degeneracy present in 
the SPT data between Yp and n^. The fact that the 
resulting Yp is higher than the BBN-consistent value re- 
flects the fact that, in the ACDM model, where Yp is re- 
quired to be BBN-consistent, SPT prefers a lower value 
of Us than does WMAP7. Since Yp has none of the late- 
time effects of neutrino species, adding BAO and Hq 
data hardly change the Yp constraint. The constraint 
from CMB+BAO+TJo is Yp = 0.305 ± 0.024. However, 
slightly more than half the parameter space allowed by 
the CMB+BAO+i/o data is ruled out by the protosolar 
upper limit. 

The constraints presented above can also be com- 
pared to the BBN predicted value of Yp = 0.24774 ± 
0.00014. The ^0.06% uncertainty quoted here does not 
include the 0.2% theoretical uncert ainty on in th e 
PArthENoPE code used in CAMB ( Pisanti et aL]|2008[ ). 
The inferred Yp value for CMB+BAO+i^o is 2.4 cr above 
the BBN prediction. We explore the relationship be- 
tween the damping-tail inference of Yp and the BBN 
prediction in more detail in § |9j where we consider an 
expanded model with both Yp and Nes free. It is worth 
repeating that the feature of the CMB data that drives 
these high values oi Yp - a trend of decreasing power at 
higher multipoles relative to ACDM predictions - can 
also be explained by other model extensions, such as 
dus/dlnk or TVeff, both discussed above. 




Fig. 14. — The CMB data place independent constraints on the primordial helium abundance, Yp. Left panel: The contours show 
the two-dimensional 68% and 95% confidence intervals for Yp and for three datasets: SPT, WMAP7, and SPT+ WMAP7. The BBN 
prediction for 3 neutrino species is marked by the horizontal black dashed line at Yp = 0.2478. The change in the damping scale due to 
varying Yp can be mimicked by changing rig for the SPT data. WMAP7 provides an independent measurement of ris, greatly enhancing 
the SPT+ WMAP7 measurement of the helium abundance. Right panel: The marginalized one-dimensional posteriors for Yp in different 
model and data combinations. We show constraints from the SPT (blue dot-dashed line), WMAP7 (blue dashed line), and SPT-(- WMAP7 
(black solid line). The BBN prediction for 3 neutrino species is marked by the vertical black dashed line at Yp = 0.2478. The orange solid 
line shows the marginally wider BBN prediction from the ACDM-fAfcff MCMC chain. The grey area is ruled out by the 95% upper limit 
from the protosolar helium abundance. The preference of the CMB for higher helium abundances can be reduced by introducing other 
extensions that affect small scale power, as is demonstrated in § |9.1| 



8.3. Early dark energy 

The dark energy density could be considerable in the 
early universe if the equation of state parameter w was 
much larger in the past, a situation commonly referred to 
as early dark energy. Early dark energy models lead to 
much larger signatures in the CMB anisotrop y than tra- 



ditional dark energy models; see|Reichardt et al. (2012a I 



for a discussion of the effects of early dark energy. One 
reason these models are interesting is that they can be 
constructed with attractor solutions that reduce the ne- 



cessity to fine-tune the initial conditions ( Wetterich 1988 
Ratra fc Peebles|1988 |. For instance in tracer" models 



the dark energy equation of state is required to be equal 
to the background equation of state - e.g., during radia- 
tion domination, w = 1/3. The advantages are that the 
DE density can stay much closer to the energy density 
of the dominant component and that the DE evolution 
is independent of initial conditions. 

Instead of considering a specific model, we choose to 
constrain early dark ene rgy in the more general tracer 
model parametrization by Doran & Robbers ( |2006 ) . This 
parametrization introduces two new parameters m addi- 
tion to the ACDM set, the dark energy equation of state 
at z = 0, Wq, and the dark energy density relative to the 
critical density at early times, ilc which is taken to be 
constant at sufficiently high redshift (z > 10). The dark 
energy density and equation of state are then given by 



w{a) - 



dc ' m 



(ilnf2dc(a) 
din a 



3(a -I- aeq) 



(15) 



(16) 



and fll 



Here a^q is the scale factor at matter-radiation equality, 
(Q'l-^) is the dark energy (matter) density relative 
to critical density at z = 0. The dark energy equation of 
state today is wq. This parametrization assumes spatial 
flatness so that -I- = 1. Since we force ^^de(a) to 
be constant at high redshift, the dark energy equation of 
state mimics that of the dominant component at early 
times, thus behaving like a tracer model. Later, during 
matter domination at z < 10, the equation of state tran- 
sitions towards its current value, wq, so it can account 
for cosmic acceleration. 

To consistently describe the perturbations, we are mo- 
tivated by quintessence models to treat dark energy as 
a perfect fluid wi th a sou nd speed, c^, equal to the 
speed of light (see Hu||1998 ). This choice, together with 
the parametrization for the background evolution, com- 
pletely specifies the behavior of dark energy. We require 
Wq > —1 and thus do not en tertain the possib ility of 
"phantom crossing" (see, e.g., Fang et al. 20081. This 
restriction allows us to avoid pathologies in perturbation 
evolution and to stay in the quintessence regime. 

The small-scale CMB temperature anisotropy power 
measurement from the SPT bandpowers improves the 
constraints on the early dark energy density over 
WMAP7 alone by a factor of 3.5; the 95% upper limit on 
r^e is reduced from 0.052 for WMAP7-only to 0.013 for 
WMAP7-I-SPT. This is a 38% improvement on the up- 



per hmit of < 0.0 18 reported for WMAP7-hKll (Re 



ichardt et al.||2012a[ ). Adding low-redshift geometrical 
measurements does not help constrain early dark energy, 
although, these data have a dramatic effect on the quality 
of the constraints on the late-time dark energy density 
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Fig. 15. — Changing N^g introduces phase shifts in the acoustic 
oscillations, however, constraints on these phase shifts are weak- 
ened by a degeneracy between A'^^ff and dg. Here, we show the 
marginalized two-dimensional posteriors for A^eff a^nd from the 
CMB data for the ACDM-|-Afcfi-l-yp model. The black curve shows 
the correlation between the two parameters due to the phase shifts 
in the a coustic oscillations. This curve is defined using the Bashin- 
|sky fc S cljak (2004j) prescription for the phase shifts, normalized 
to the best-fit 6»s = 1.0429 X lO'^ at N^g = 3.046. 



and equation of state. The upper limit is essentially un- 
changed at < 0.014 for WMAP7+SPT+BA0-f SNe. 
The < 0.013 bound from WMAP+SPT is the best 
published constraint from the Cosmic Microwave Back- 
ground. 

9. TWO-PARAMETER EXTENSIONS 

We will now consider two 2-parameter extensions to 
ACDM: Ncs+Yp and N^ff+Y, m^, . Of the many pos- 
sible 2-parameter combinations, we limit our discussion 
to these two extensions because they are physically well- 
motivated. The first case is an interesting consistency 
test of BEN, while the second case is a natural space to 
consider for sterile neutrinos. 



9.1. iVoff and Yp 

Here we examine simultaneous constraints on the effec- 
tive number of neutrino species and the primordial abun- 
dance of ^He: ACDM-l-A'cff-l-1^. This two-parameter ex- 
tension to ACDM is moderately favored by the data with 
a Ax^ = 5.5 (see Table [2|, though it does little to ease 
the tension between datasets beyond the one-parameter 
ACDM-I-A^eff model. 

We first discuss the physical mechanisms that lea d to 
joint constraints on A'^ff and Yp. As noted earlier in 



and |8.2[ both the number of neutrinos and helium abun- 
dance primarily impact the damping scale. However, the 
parameters are not fully degenerate due to se veral addi- 
tional effects of varying A'eff as discussed in Hou et al. 
(2011). Here we highlight in particular the role of shifts 
in the acoustic peak loca tions induced by phase shi fts in 
the acoustic oscillations ( Bashinsky fc: Seljak||2004 ). 

Acoustic oscillations in a constant gravitational poten- 
tial will have the form cos{krs{r]) + (p), where k is the 
wavenumber and the phase = as a consequence of 



the initial conditions in inflation models. Changes in the 
gravitational potentials shift the phase away from zero. 
Unlike photons, neutrinos have a long free-streaming dis- 
tance in the early plasma. They thus alter the evolution 
of the gravitational potential and therefore the resulting 
phase shifts. For modes that entered the horizon in the 
radiation-dom inated era {I > leg — 434), Bashinsky & 
Seljak (2004) find that the change to the phase shitt due 
to neutrino free streaming is A(p — 0.19Tr pi,/{piy + p~^). 
Here and pi, are the energy densities of the photons 
and neutrinos respectively. This phase shift changes the 
positions of the acoustic peaks in the CMB data. 

Constraints on the phase shifts are weakened by a par- 
tial degeneracy between the angular size of the sound 
horizo n an d a phase shift. This point is illustrated by 
Figure [15] which shows the 68% and 95% confidence in- 
tervals on A"off and 9s from the CMB. Changing the phase 
shift by 5tf and the sound horizon by 59s moves the m*^ 
acoustic peak by 5lm = {Sf — lmS9s)/0s- We can thus 
hold fixed the position of any single acoustic peak by 
varying 9s to accommodate the phase shift. The black 
curve in the figure is the locus of points that does ex- 
actly this, although not for a particular peak but for a 
particular / value. We chose / = 1666 because it is ap- 
proximately where the signal-to-noise is maximized in 
the SPT bandpowers. The major axis of the constraint 
ellipse falls along the black line, confirming that the cal- 
culated phase shifts explain the observed correlation be- 
tween A"cff and 9s. The degeneracy between A^off and 9s 
is broken (the contours close) in part because the com- 
pensating value of 9s is Z-dependent; there is no one value 
of 9s that preserves all the peak locations. Other effects 
play a role in the degeneracy breaking as well, including 
the early ISW effect and effects due to the high values 
of the baryon fraction that occur at low A"off (Hou et al. 
20TT|). 



The Ncff-Yp contours from the CMB dataset are shown 
in Figure [I6] The contours form an ellipse, with the mi- 
nor axis (the better constrained direction) oriented in the 
direction of varying 9d/9s. Relative to WMAP7+K11, 
the new SPT bandpowers have the most impact on the 
high Aoff' region. This improved constraint on N^g is due 
to improved constraints on the neutrino-induced phase 
shifts. 

The inferred values of Ncff and Yp for the CMB data 
are: 



Ip^ 0.314 ±0.033, 
Afeff = 2.60 ±0.67. 



(17) 



When the BAO and Hq observations are added, the me- 
dian values of these parameters shift slightly towards 
the standard model expectations, and the uncertainties 
tighten: 

(18) 



Ip^ 0.294 ±0.030, 



A^, 



off' - 



: 3.32 ±0.45. 



We compare these CMB-derived results with the pre- 
dictions from BBN theory (see Equation 12 ) in Figure 16 
The dashed curve in Figure [16] denotes iKe BBN theory 
predictions. The data are consistent with BBN theory, 
although they display a mild preference for higher-than- 
predicted values of Yp. The region above the 2a upper 
limit from the initial solar helium abundance (Serenelli 
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Fig. 16. — This plot contrasts the CMB-derived two-dimensional 
likelihood contours for N^fi and Yp with the predictions of BBN 
and other Yp measurements. The black solid contours are the 68% 
and 95% confidence intervals for SPT+WMAP7. The scattered 
points are samples from the Markov chain with the color encoding 
the value of 9^/ 9s. The dashed black curve marks the standard 
BBN rela tion ship between A^^ff and Yp. Similar to the right panel 
of Figure [14] the gray area is ruled out by the 95% upper limit of 
the measurement of solar initial helium abundance by Scr cnelli ff| 
[Basu (2010). The Yp measurement from spectral observations of 
extr agalact ic low-meta l licity ionized gas clouds by |Peimbert et al.| 
([2007| and |Aver et aT] | |2012[ | are shown by the light blue and red 
bands respectively. 



&: Basu 2010 ) is gre yed out. We also plot the Peim- 
bert et al.l ('2UU7) and I Aver et all (|2012D estimates ot 



based on spectroscopic ob serv ations ot low-metallicity ex- 
tragalactic clouds (see § 8.2) as horizontal bands. The 
CMB-|-BAO-|-i?o inference tor the helium abunda nce is 
consistent w ith the spectroscopic observations of |Aver| 



et al.| ( |2012[ ) at 1.3 cr 

In summary, the CMB data can constrain Yp and 
Ncff simultaneously by measuring the damping scale and 
the location of the acoustic peaks. The combination of 
CMB+BAO-|-ifo prefers a value of N^s which is con- 
sistent with three neutrinos species and a value of Yp 
slightly above the BBN prediction. 

9.2. Ncff and ^ 

We finally focus on a two-parameter extension, 
ACDM-|-^ TOiy-l-A^off , in which both the mass and num- 
ber of neutrino species are freed. This is a well- 
motivated, natural model expansion of the individual 
J2 and NeS extensions discussed earlier since (1) neu- 
trino oscillation ex periments have establishe d that neu- 
trinos are massive ( Hamish Robertso n|2009 1 and (2) the 
event excess in the 



MiniBooNE search tor v, 



oscil- 
lations (lAguilar-Arcvalo et al. 2010) and the reactor an- 



tineutrino anomalies (^Mention ct al.^2011^ suggest one or 
more massive (~ 1 eV) sterile neutrmo species. We will 
continue to assume equal masses in all neutrino species. 
For treatment s of (3-1-1) or (3-1-2 ) m odels in the litera- 
ture, see e.g. Kopp et al. ( 2011[ ) or Giunti & Laveder 



Cosmological observations ca n rule out large ar 



eas of the Nes-^rrii, plane. [Benson et al.| ([20111 
investigated constraints from WMAF7, Kll, BAU, 
Hq, and SPT galaxy clusters, and found Ncs = 
3.91 ± 0.42 and J2mt, = (0.34 ± 0.17) eV. A 



used WMAP7-hKll-hWiggIeZ+7TlTI+SNrSB"^to mea 



sure iVoff — 3.58 



+0.15 
'0.16 



CL) and 



< 



0.60 eV (95% CL). The tightness of the constraints on 
Ncs from this study are driven by the inclusion of the 
H(z) inference from measuring the age-redshift relation 
of passively evolving galaxies out to z ~ 1.75 (Moresco 
et al.[[2012|), which depends heavily on modeling stel- 



lar evolution. The extended arm out to higher redshifts 
provides greater constraining power on the expansion his- 
tory, placing constraints on TVeff! conversely, this infor- 
mation has little impact on the neutrino mass constraint. 

As discussed in § [4.2[ this two-parameter exten- 
sion significantly improves the quality of fit to the 
CMB+BAO+i/o data by = 7.9, equivalent to a 2.3 cr 
Gaussian preference. This model also increases the con- 
sistency between the CMB, BAO, and Hq datasets, as il- 
lustrated by Figure 17 Both panels show the rs/Dv{z = 
0.57)-7Jo plane and~are identical except for the color- 
coding, which encodes A'off in the left panel and ^ 
in the right panel. Like the ACDM-I- A'off model shown in 
Figure [31 this 2-parameters extension allows the contours 
from the three datasets to overlap comfortably. 

iVgfj and ^ nil, are individually constrained by sepa- 
rate features in the CMB as described in § [6| and § [8.1| 
- the early ISW effect for ^ m^^ the damping scale and 
the position of the acoustic peaks for A'off. We can use 
Figure [T7| to understand phenomenologically the relative 
interplay between A^eg, Hq, and the BAO fea- 

tures, given a CMB prior. In the left panel, A'off in- 
creases as Hq increases and rg/Dy decreases. The de- 
tailed direction means that the Hq measurement is much 
more constraining than BAO for A'oft, si mila r to the con- 
straint on individual A'off discussed in § [8.1| In the right 
panel, the BAO data are critical for the neutrino mass 
measurement, whereas varying Hq at fixed rg/Dy (hori- 
zontal lines in the figure) has almost no effect on ^m^. 
This can also be seen by the significant tightening of the 
inferred mass unc erta inty between the gre en and blue 
contours in Figure [18] We see in Figure [Ts] that there is 
a positive correlation between A'off and 2J when Hq 
is one of the datasets. This correlation emerges because 
the only way to increase A'efi in this model space while 
fixing Zoq, Os and Hq is to increase '^rrii, (and decrease 

The CMB constraints on A'efi and ^ rrii, are nearly 



independent as shown by the red contours of Figure 18 
As would be expected for independent parameters, the 
quality of the constraints are nearly unchanged from 
the single-parameter extensions discussed earlier. The 
A'off constraint is A'off — 3.40 ± 0.48 with massive neu- 
trinos compared to Aoff ~ 3.62 ± 0.48 in the ear- 
lier massless case. In both the ACDM-|-^ nii, and 
ACDM-^Aoff-^X) "ii/ models, the CMB sets an upper 
limit on the neutrino masses of < 1.6 eV (95% 

CL). This independence is broken, and the constraints 
significantly tightened, once low-redshift observations 
are added. The green contours show the results of 
including Hq, while the blue contours show the con- 
straints from CMB+BAO+ffo- For the combination of 
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Fig. 17. — The Hq — rs/Dy plane for the two-parameter extension ACDM-|-7V(,£f+ m^^ with the color of scattered points coding the 
values of N^g (left panel) and X]m^ (right panel). The colored contours and the gray filled regions are as described in Figure [s] 

This is a greater than 3 a preference for ^ mi, > and 
a 2.2(7 preference for A'off > 3.046. 
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Fig. 18. — This figure demonstrates the impact of each combina- 
tion of datasets on the constraints on ^2 '^"^ ^^^d ^cff . The shaded 
contours are the 68% and 95% confidence intervals for the fol- 
lowing data combinations: SFT+WMAP7 (CMB; red). CMB+Ho 
(green), CMB-f Ho+BAO (blue), CMB-f //o-f BAO-(-SPTcl (pur- 
ple). The combined data are in >2a tension with the ACDM 
assumption of three massless neutrino species. 



CMB+BAO+Ho we find: 

^ m^ = (0.48 ± 0.21) eV, (19) 

Aoff = 3.89 ±0.37. 

As should be expected from § [6j the mass constraint can 
be improved by adding a tracer of structure growth. For 
CMB+BAO-fHo+SPTcL (the purple contours, we find 
the tightest constraint of 



^ = (0.51 ± 0.15) eV, 
^eff = 3.86 ± 0.37. 



(20) 



10. CONCLUSIONS 



The SPT bandpowers as presented in |Story et al.| 
(20121 are currently the best measurements of the CMB 
damping tail from the third to ninth acoustic peaks 
(650 < / < 3000). The SPT bandpowers greatly improve 
measurements of three CMB observables: the damp- 
ing scale due to photon diffusion, the locations of the 
acoustic peaks, and the amplitude of the lensing poten- 
tial. The combined CMB data show a slight, reddening, 
scale-dependent tilt unreproducible in the ACDM model. 
Therefore, we find that while the SPT bandpowers are 
well-fit by a ACDM cosmology, there are tantalizing hints 
for extensions to this model. 

With the aid of Table [5] we now summarize our most 
intriguing findings. The extension that the CMB data 
most prefer (at 2.4 a) is the running of the spectral in- 
dex, dris/dlnk. This extension improves the fit to the 
CMB data, reducing the minimum by Ax^ — 4.9 
from that achieved with ACDM. We find, for the CMB 
alone, dng/dlnk = —0.024 ± 0.011. Negative running 
of the spectral index reconciles the slightly redder spec- 
tral index preferred by the SPT data with the slightly 
bluer tilt preferred by the WMAP7 data; the same value 
of running is preferred by the WMAP7 data alone, al- 
though only at ^ la. The combination of CMB, BAO, 
and Hq data prefers more negative running, dus/dhik — 
—0.028 ± 0.010. Deviations from scale invariance of this 
magnitude are not expected in standard inflationary sce- 
narios. 

Increasing the primordial helium abundance increases 
the damping scale which has a similar impact as neg- 
ative running on the CMB power spectrum. Freeing 
the helium abundance reduces the minimum of the 
combined CMB fit by Ax^ = 4.4, resulting in 1^ = 
0.300 ± 0.025. However, this value of Yp exceeds the 
protosolar 95% confidence upper limit of Yp < 0.294. A 
summary of which extensions the CMB prefer can be 
found in Table [H 

We also find that some extensions are more motivated 
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TABLE 5 
Hints for extensions to ACDM 



Model extension 


Preferred by the 
CMB data 

(see 1=* row of Table [Sjl 


Alleviates tension 
when combining 
BAO and Ho data 


Statistical error on extension 
parameter reduced by combining 
BAO and Hq data 




Marginal 

Ax^ = 2.4, 1 dof 

Figure |9] 


No 

Figure [3] 


Yes 

Figure |9| 


\0+Ho/°'CMB = 0.37 


dus/dlnk 


Yes 

Ax^ = 4.9, 1 dof 
Figures [Tl| and [Ti] 


No 

Figure [i] 


No 

""CMB-IjEj 
Figure |l] 


JO+Ho/°'CMB = 0.91 




No 

Ax^ = 1.1, 1 dof 
Figure |l3| 


Yes 

Figure [s] 


Marginal 

""CMB-tEi 
Figure [13 


iO+Ho/^'CMB = 0.73 

] 




No 

Ax^ = 2.6, 2 dof 
Figure [lil 


Yes 

Figure [Tt] 


Yes 

o'Cmb+bao+Hq/o'Cmb = 0.43 (Y^m,^) 
fCMB-taAO+iTa/o'CMB = 0.77 (A^off) 
Figure [l8| 



This tabic summarizes the evidence for physics beyond the ACDM model. We state our conclusions and 
relevant figures for each entry. The second column addresses whether the CMB alone shows a preference 
for the extension. The third column addresses whether the extension reduces the tension between the 
datasets. The Ax^s in Table |2| quantify the preference and reduction in the tension. The final column 
shows if the extension is favorea because the statistical errors are significantly reduced by adding datasets. 



than others with respect to how they reduce the mild ten- 
sion between the CMB, Hq, and BAO datasets in the con- 
text of the ACDM model. The ability of various exten- 
sions to reduce this tension is shown in the third column 
of Table 5. Of the extensions that we consider, the effec- 
tive number of neutrino species, N^s, is the only model 
extension that significantly reduces the tension between 
the CMB, BAO, and Hq datasets. It does so both on its 
own, and when ^ m,y is simultaneously free. The tension 
is reduced because iVefi brings the model predictions for 
CMB and into better agreement with data without 
much change in the prediction for the BAO observable 
Vs/Dv- Thus, this extension can bring the region of the 
CMB probability posterior that is more consistent with 
the measured Vs/Dy toward simultaneous compatibility 
with the Hq measurement. Combining the CMB with 
BAO and Hq data, we find N^e = 3.71 ± 0.35. Simul- 
taneously freeing ^ m^, leads to further improvement in 
the fit as higher neutrino masses bring the CMB into bet- 
ter consistency with rs/Dy, while increasing Ncs leads 
to consistency with Hq. 

The impact of combining low redshift measurements 
with CMB data is summarized for selected extensions 
in the last column of Table 5. We find that the CMB, 
CMB+BAO, and CMB+SPTcl aU at least weakly prefer 
a non-zero neutrino mass (with only Hq preferring zero 
mass). The combined CMB-hBAO-h^^o+SPTcL dataset 
constrains ^ = (0.31 ± 0.11) eV. The preference for 
massive neutrinos remains at > 2 cr as long as we keep 
at least one of the two most precise BAO measurements; 
however, it disappears if we do not include any BAO 
data. 

We find improved constraints on two additional model 
extensions, early dark energy and non-zero curvature. 
These extensions do not address the scale-dependent tilt 
present in the data or lead to a reduction of the tension 
between the CMB and low-redshift data. 

In summary, we explore six theoretically motivated ex- 



tensions to the ACDM model, and find that the combi- 
nation of the SPT data with other cosmological mea- 
surements yields a 2-3 a preference for some of these 
extensions. Favored extensions are running of the pri- 
mordial spectral index, a sum of neutrino masses of or- 
der 0.3 eV, or greater than 3 effective neutrino species. 
The evidence for any extension to the standard ACDM 
model is currently weak. However, a significant detec- 
tion of any of these extensions in future data sets such 
as Planck and high-resolution CMB polarization exper- 
iments would dramatically impact our understanding of 
particle physics and cosmology. 
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